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Abstract 

In recent years, time-dependent quantum mechanical approach has been used exten- 
sively to investigate various elementary chemical processes such as photo dissociation, pre- 
dissociation, reactive scattering, molecule-surface scattering etc., In the present study we 
have investigated (i) photo-predissociation of Nal using femtosecond laser pulse, (ii) chan- 
nel control in chemical reactions using HI photodissociation as a case study and (iii) 
photodissociation and predissociation processes in OH radical. 

In chapter 1 of the thesis, the necessary background material related to the TDQM 
approach to photodissociation and predissociation processes have been reviewed. We have 
also presented a short literature review emphasising mainly on the recent studies. 

In chapter 2, theoretical and computational methodologies related to the TDQM ap- 
proach to photo dissociation and predissociation processes have been presented. This chap- 
ter has been divided into three sections. In section 1, we have described the existing com- 
putational methodologies for implementing the TDQM approach. Section 2 of this chapter 
is devoted mainly for the interaction of light with molecules. We have presented a semi- 
classical and a quanttun mechanical description of the photoabsorption process. We have 
also described, how the photodissociation cross section is computed in the time-domaiji. 
Also we have presented a critical analysis of the nature of the wavefunction prepared by 
a laser pulse on the excited electronic state. In section 3, we have described the computa- 
tional methodologies for investigating multieletronic state problems. Also the golden rule 
treatment for computing the predissociation linewidth is included in this chapter. 

In chapter 3, we have presented the results of the TDQM study of the photo- pre- 
dissociation of Nal using ultrashort laser pulse. We have investigated the nature of the 
wave packet prepared by a femtosecond laser pulse, whose excitation frequency is chosen 
near, very near and far from resonance excitation frequency. Since the excitation fre- 
quencies are different, the WPs prepared on the excited electronic state also have different 
characteristics. The norm of the WP on the excited electronic state during the pulse 
preparation step has been analysed. After the pulse preparation step is over, the WP is 
time-evolved for about 3 ps on the covalent excited electronic state which is coupled to 
the ground electronic state. Due to the diabatic coupling between the ionic groimd state 
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and the covalent excited state, the time-evolving WP bifurcates near the crossing point 
and a small fraction of it (as determined by the Landau-Zener curve crossing probability) 
leaks through the crossing to form neutral products Na and I and the remaining part of 
the WP oscillates between the ionic and the covalent turning point. Whenever the WP 
reaches near the crossing point a small fraction of it leaks through crossing. We have 
estimated the oscillation periods corresponding to photoexcitation wavelengths 300, 321 
and 344 nm and our results are in good agreement with the experimental results of Zewail 
and coworkers. However, when the excitation wavelengths are far from resonance, i.e., Aex 
= 360 and 390 nm, the WP spreads immediately ( in less than 0.5 ps ) and shows several 
interference patterns. Therefore, the oscillation period estimated by us for \ex = 360 and 
390 nm are not in agreement with experimental results. In order to simulate them, a two 
photon calculation is required as the experimental results are strongly dependent on the 
probe wavelength. We have observed that the transient effect and the energy of the WP 
are highly responsible for the unusual characteristics of the WP evolution. 

In chapter 4 of the thesis, we have presented a detailed analysis of the photodissociation 
of HI in its first absorption band. HI photo dissociation takes place around 250nm resulting 
in H -f I and H -h I* formation. In this energy range, the following potential energy curves 
participate in the photo dissociation: ^IIi, ^Hi, and ^Hq. Among these, the former two 
correlate with the H -h I channel and the later two correlate with the H -t- 1* channel. Except 
for the ^Ho state, all other states are non-adiabatically coupled to each other. We have 
computed the photodissociation cross section(cr) corresponding to the photoexcitation from 
the ground vibrational level of the ground electronic state to all of the above mentioned 
excited states. We have made use of the PE curves, non-adiabatic coupling potentials and 
the coordinate dependent transition dipole moments reported earlier by Levy and Shapiro. 
Our computed a and E values are in excellent agreement with the experimental and the 
time-independent quantum mechanical results. We have also computed the anisotropy 
parameter(/3) for the H -I- I* channel and our /? values as a function of excitation energy 
are in accord with the experimental and the time-independent quantal results. In addition, 
we show that the P value changes significantly with vibrational excitation of the parent 
molecule and point out vibrational excitation prior to photodissociation as a means of 
channel control. The results are interpreted in terms of the influence of the non-adiabatic 
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coupling on the time-evolution of the wavefunction. 

In chapter 5, results obtained for the photodissociation and pre dissociation processes 
in OH radical using the TDQM method ‘ have been described. Since OH radical is one 
of the important constituents in the astrophysical environment, understanding various 
dissociation pathways in OH becomes particularly essential. Photo dissociation in OH 
takes place due to electronic transition from the ground electronic state (X^H) to the 
excited electronic states, I'A, and 2^TL. The PE curves and the coordinate 

dependent transition dipole moments are taken from the ab initio results of Dalgarno and 
coworkers. Our computed photodissociation cross sections are in good agreement with the 
time-independent quantal results and a factor of two less than the experimental results. 
Predissociation in OH takes place from its first absorption band (X2n-A2E+) due to spin- 
orbit coupling between the boxmd state A^E'*' and the repulsive states '^E , l^E and “^H. 
Using the recent ab initio results of PE curves and the spin-orbit matrix elements reported 
by Yarknoy we have obtained predissociation linewidths corresponding to predissociation 
from v’=l-7 of the A^E^ state. The results are analysed in terms of the crossings between 

the PE curves. 

A summary of our findings are given in chapter 6. 
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Chapter 1 


Introduction 

Photon-induced processes have drawn considerable attention in recent years[l-5]. One 
of the most important photon-induced processes which offers a wide variety of applications 
is photodissociation{&^. In chemical dynamics, for instance, photofragments produced by 
photodissociation are used in investigating bimolecular exchange reactions in molecular 
beam experiments[7-10]. Photodissociation processes have technological applications as 
well[ll, 12]. In addition, a careful study of photo dissociation provides an opportunity 
to verify the existing dynamical theories. Besides these applications, photodissociation 
is one of the fundamental processes that occur in earth’s atmosphere and astrophysical 
environment [13-15]. Therefore, understanding the dynamics of photofragmentation be- 
comes particularly important. 

Photodissociation processes usually involve electronic transition (due to photoabsorp- 
tion) from the ground electronic state to an excited repulsive electronic state. As a result, 
the molecule fragments into its constituents [16]. For a diatomic molecule AB, this process 
can be represented as 

where v and j refer to the vibrational-rotational state of the parent molecule, a; is the 
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frequency of the photolysing radiation. Unlike a full collision process ( e.g., a bimolecular 
reaction), the short lived intermediate ( or the transition state, [AB]*'^) is produced via 
photoabsorption. Thus we are mainly interested in the second half of a full collision event. 
Therefore, photodissociation process is often referred to as a half-collision process[17]. 
For such a half-collision process, one can experimentally measure scalar quantities such 
as photodissociation cross section and energy distributions of the photofragments and the 
vector quantities such as angular distribution of the products, anisotropy parameter and 
angular momentum correlations[18-25]. 

Sophisticated experimental techniques developed in recent years enable us to obtain 
these quantum state resolved quantities precisely. For instance, one can prepare the 
reagent in specific vibrational-rotational quantum state by techniques like stimulated emis- 
sion piunping or infra-red laser excitation prior to photodissociation[22-25]. Similarly 
photofragment velocities can be measured by time-of-flight mass spectrometry and the 
internal state distributions of the products can be obtained by laser-induced-fluorescence 
(LIF) technique or resonance enhanced multiphoton ionisation (REMPI) technique[21-25]. 
The photofragment imaging technique developed recently by Chandler and coworkers [26, 
27] enable an aw;curate determination of photofragment angular distribution. Zewail and 
coworkers [29-32] have demonstrated that ultrashort (femtosecond) laser pulses could be 
used to map out the molecular motion on the excited electronic state in real time. Thus 
the experimental methods have reached the ultimate time resolution[32] . 

Theoretical developments have kept pace with experiments. Within the framework 
of quantinn mechanics, one can obtain the quantum state resolved cross section either 
by time-independent quantum mechanical approach (TIQM) [35-42] or by time-dependent 
(TDQM) [43-65], In the former approach, the photodissociation cross section is computed 
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from Franck-Condon factors. Therefore, one has to know the initial and the final scatter- 
ing state eigenfunctions in order to compute the cross section. However, evaluation of the 
scattering state eigenfunction is computationally more demanding especially when there 
are more than one product channel open. Such problems are circumvented in the TDQM 
approach [66]. The time-domain formulation of spectroscopy was pioneered by Heller and 
coworkers [43,67-74]. The computational simplicity of this approach and the physical in- 
sight offered it are primarily responsible for the success of the TDQM method[75, 76]. 

Heller and coworkers [75-77] employed the Gaussian wave packet (GWP) method ( a 
semiclassical treatment) for following the time-evolution of the system. The GWP method 
assumes that the potential-energy-surface (PES), upon which the dynamics occurs, is lo- 
cally quadratic and that the process of interest takes place in short time such that the 
Gaussian wave packet remains Gaussian during the time-evolution [75-81]. Since direct 
photodissociation and resonance Raman scattering processes occur usually in less than 
100 fs, Heller and coworkers could successfully implement the GWP method for a variety 
of systems[43, 67-74]. Until 1983 there was not much progress in using the TDQM ap- 
proach as there were not efficient numerical algorithms and computer memory and speed 
were limited. However, the situation has changed now. Computers with large memory 
and high speed have become common. The fast Fourier transform(FFT) method [82-86] 
is being used extensively to obtain the spatial derivatives of the wavefunction and it is 
found to be accurate and much faster than the methods employed earlier [83, 85]. Several 
time-evolution methods have been developed in the last decade or so. The numerically ex- 
act WP propagation method (often referred as the grid method) offers several advantages 
over TIQM approach [87-90]. In addition to providing the computational simplicity and 
the physical insight, a single wave packet propagation is sufficient to obtain photodissoci- 
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ation cross section over entire photoexcitation energies. Therefore, the numerically exact 
WP propagation method has been used rather widely in recent years for studying several 
elementary chemical processes[87, 88]. For an extensive review of them and their relative 
merits, the reader is referred to several reviews that have appeared[87, 88, 91-97]. Several 
reviews [87, 88, 93] and monographs[95-97] on the applications of the TDQM approach 
have emerged in recent years. 

One of the major difficulties with the photoabsorption and related processes is that 
there could be more than one electronic state participating in the dynamics[98-113]. There- 
fore, one has to go beyond the Born-Oppenheimer approximation and use either adiabatic 
or diabatic representation to investigate the dynamics[98, 115, 116]. In principle, one can 
extend the existing time-evolution methods for investigating the dynamics which takes 
place on multiple electronic states[117-120]. However, Coalson and Kinsey [78, 121, 122] 
introduced the WP perturbation theory approach and Alvarellos and Metiu[123, 124] in- 
troduced a computational algorithm by making use of both adiabatic and diabatic rep- 
resentations to study especially the coupled electronic state problems. Thus, using the 
quantum mechanical approach, quantum state resolved cross sections (partial photodis- 
sociation cross sections [36, 43, 54, 65, 90, 125]), branching ratio between the difiFerent 
products when there is more than one product channel [79, 126] and angular distribution 
of the products can be determined accurately [20, 35, 127-130]. Apart from these devel- 
opments, theoretical and computational studies have emerged recently to investigate the 
nature and the dynamics of the WP prepared by an ultrashort laser pulse[131-137]. Sev- 
eral theoretical studies have also been carried out to simulate the experimental studies of 
Zewail and coworkers[132, 138-147]. 

In addition to computing cross sections and understanding the dynamics, an investi- 
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gation of photodissociation also provides clue to means to control the product formation 
ratio[148, 154, 155]. Recently, various channel control methods which are mainly based on 
photoexcitation, have been developed to influence the product formation ratio [149-1 66]. 

Over the years, photodissociation in several diatomic and polyatomic molecules have 
been studied using the TDQM approach. A review of literature covering up to 1992 have 
appeared elsewhere[3]. Therefore, we only cite the most recent studies in table 1.1. 

Another important photon-induced process which deserves attention is the predis- 
sociation [16, 116]. In predissociation process the photoabsorption takes place from a 
bound (ground) electronic state to a bound excited electronic state (bound — + bound transi- 
tion). Before the excited molecule emits the light and returns back to the ground electronic 
state, it dissociates because of coupling between the bound excited electronic state and a 
repulsive excited electronic state. Depending upon the perturbation present in the system, 
predissociation is classified as electronic, vibrational, rotational, spin-orbit and spin-spin 
types[116]. Except for rotational predissociation, other types of predissociation involve 
participation from more than one electronic state. Since predissociation lifetimes are of 
the order of 10~®-10“^^ seconds, there are only limited number of TDQM studies available 
in the literatme [167-1 70]. However, vibrational predissociation in van der Waals com- 
plexes have been studied extensively in recent years as they provide mechanistic insight 
into the intramolecular vibrational redistribution process[171, 172]. Some of the recent 
studies on predissociation are also listed in table 1.1. 

In the present investigation, we have undertaken a TDQM study of photodissociation 
and predissociation processes in some diatomic molecules. Since diatomic molecules are 
relatively simpler systems, a detailed insight into the dynamics of fragmentation can be 
obtained. As we have mentioned earlier, nonadiabatic coupUng among potential-energy 
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curves are inevitable in the excited state dynamics. Therefore, we have paid special atten- 
tion to this phenomenon. We have investigated the nature and the dynamics of the WP 
prepared by an ultrashort laser pulse by taking Nal predissociation as a case study. We 
have analysed in detail the photodissociation of HI in its first absorption band. We have 
also investigated the effect of vibrational excitation of HI prior to photodissociation on 
the product formation ratio (i.e., I* /I branching ratio). Also we have studied photodisso- 
ciation and predissociation processes in OH radical. The thesis has been divided into six 
chapters. In chapter 2, theoretical and computational methodologies for implementing the 
TDQM method to study photodissociation and predissociation processes are described. In 
chapter 3, Nal predissociation using femtosecond laser pulse is presented. In chapter 4, 
we present the results of our studies on HI photodissociation. In chapter 5, OH radical 
photodissociation and predissociation processes are described. Summary and conclusions 
are given in chapter 6. 
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Table 1 . 1 : Summary of recent TDQM studies on photodissociation/predissociation pro- 
cesses 


SI. No. 

Year 

System 

Time- 

evolution 

Laplacian 

Reference 

1 

1993 

CH 3 I 

SOD 

FFT 

[1731 

2 

1993 

CH 3 ONO 

Chebyschev 

FFT 

[174] 

3 

1993 

HBr/LiF(001) 

Chebyschev 

FFT 

[175] 

4 

1993 

HCO 

Chebyschev 

FFT 

[176] 

5 

1993 

Na 2 

SOP 

FFT 

[177] 

6 

1993 

IBr 

Lanczos 

FFT 

1119] 

7 

1993 

HOCl 

Chebyschev 

FFT 

[178] 

8 

1993 

CH 2 

Lanczos 

FFT 

[179] 

9 

1993 

HF 

SOP 

FFT 

[180] 

10 

1993 

HI 

SOD 

FFT 

[ 120 ] 

11 

1993 


SOP 

FFT 

[181] 

12 

1994 

HBr/LiF(001) 

Chebyschev 

FFT/DVR 

[182] 

13 

1994 

CI 2 

SOP 

FFT 

[183] 

14 

1994 

Ar...HCl 

Chebyschev 

FFT/DVR 

[184] 

15 

1994 

OH 

SOD 

FFT 

[170] 

16 

1994 

O 3 

Lanczos 

FFT 

[185] 




Chapter 2 


Methodology 

This chapter has been divided into three major sections. In the first section we address 
the question, how the dynamical evolution of a system in the absence of any applied field 
can be followed quantum mechanically. The second section is devoted to a discussion of 
the interaction of radiation with matter and the resulting processes. In the third section, 
we describe how the dynamical evplution of the system is followed on nonadiabatically 
coupled electronic states. 

2.1 Time-dependent quantum mechanical approach 

The dynamical evolution of any quantum system can be followed by solving the time- 
dependent Schodinger equation (TDSE) 

“ ■&(»)*(*) ( 21 ) 

where H{t) is the Hamiltonian operator of the system. Since this is a first order differential 
equation, its solution will be of the form 

^(t) = Uit^toMk) (2-2) 
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where U(t,to) is the time-evolution operator. It is a linear operator and it satisfies the 
initial condition 

U{to,to) = i (2.3) 

For a Hamiltonian which is explicitly time-independent, the formal solution of eq. (2.1) is 
given as 

^(i) = (2.4) 

By comparing eq. (2.4) with eq. (2.2), the time-evolution operator is identified as 

U{t,to) = (2.5) 

Hamiltonians with explicit time-dependence will be discussed in the later paurt of this 
section. 

In what is termed the grid method, solving eq. (2.1) involves the following steps[85, 88, 
87, 186]: (i) representation of the initial wavefunction on a finite grid, (ii) propagation of 
the wavefunction till the process of interest is over and (iii) analysing the final wavefunction 
for dynamical observables. 

For a photodissociation process, the initial wavefunction is taken as the wavefunction 
corresponding to the state "prepared "(see below). The grid points should be close enough 
to contain all the oscillations of the wavefunction during its evolution. The grid spacing 
also determines the maximiim energy that can be represented on the grid and it should 
be chosen according the energy spectrum of the Hamiltonian. The grid size (and spacing) 
chosen in the coordinate space is inversely related to the size in the momentum space. 
Therefore, if the initial wavefunction is narrow in the coordinate space, it will be spread 
out in the momentum space. These considerations are particularly important as the Hamil- 
tonian operator consists of kinetic and potential energy parts and the former is local in 
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the momentum space while the later is local in the coordinate space. Thus one takes ap- 
propriate grid spacing and optimal width for the initial wavefunction wherever possible. 
Once the initial wavefunction is chosen, the next step is its propagation in time. As we 
can anticipate from eq. (2.4) the time-propagation will require repeated operations of H 
on 'S'. 

2.2 Evaluation of 

The evaluation of consists of the evaluation of the kinetic energy part, and the 
potential energy part, V^. As we mentioned earlier in the text, the latter is local in the 
coordinate space and its evaluation amounts to a mere multiplication of ^ at each grid 
point by V. However, the kinetic energy operator is nonlocal in the coordinate space and 

its evaluation is more complicated. A similar problem would arise with Vin the momentum 

^ ^ 

space. Therefore, it would be better to solve the and W parts of the Hamiltonian in the 

representation in which each is local. The coordinate and momentum spaee wavefunctions 
are related through the Fourier transformation as 


^(r) = 

1 

\/2^ • 

r ^{k)dk 

J—oo 

(2.6) 

m = 

1 

r e-*^ ^{r)dr 

/— OO 

(2.7) 

^'(r-) = 

1 

• 

r e*^ (ik) ^k)dk 

l~oo 

(2.8) 

ijr"(r) = 

1 

roo 

/ {ikf ^k)dk 

1—00 

(2.9) 


Therefore, the evaluation of kinetic energy operator can be carried out by the following 
steps: [85, 90] 

1. the wavefunction ^{r) is (forward) Fourier transformed into momentum space 
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2. the momentum space wavefunction is multiplied by ^ at each momentum 
grid point and 

3. the momentum space wavefunction is brought back to the coordinate space by 
an inverse Fourier transformation. 

The availability of fast Fourier transform (FFT) algorithms [187] enables us to accu- 
rately evaluate 7^(r). This method requires the wavefunction to satisfy periodic boundary 
conditions and for band-limited functions it is exact. In addition, the numerical grid should 
have optimal number of equally spaced sampling points as decided by the phase-space 
volume[87]. An attractive feature of t^e FFT method is its computational scaling prop- 
erty with respect to the number of grid points N; computer time increases as ~ NlogN. 
Usually N is taken as a power of 2 as most available FFT routines reqtdre this. It is 
worthwhile to mention that the Fourier method can be parallelised easily and hence the 
computer time scaling can be reduced as ~ (logN). 

In the discrete representation, the operators corresponding to the original Hamiltonian 
are mapped onto a discrete Hilbert space. The mapping has to be carried out in such 
a way that the operators in the discrete Hilbert space obey all the quantum mechanical 
commutation relations obeyed by the corresponding physical observables in the original 
Hilbert space. Only then one Can say that a one-to-one mapping of the operators to the 
discrete space is achieved. Kosloff and Kosloff have shown[85] that this is indeed the case 
with the Fourier discretisation for a band-limited function with a finite support - that is, 
of finite extent in configuration as well as momentmn space. This once again reveals the 
versatility of the Fourier transform method over the other methods. If the wavefunction is 
not band-limited then one has to employ a semiclassical approximation for kinetic energy 
operator as has been done in the finite difference scheme. 
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2.3 Time-evolution 

Having represented the wavefunction and the Hamiltonian on a numerical grid, the next 
step in the TDQM approach is to propagate the wavefunction in time. This heis to be 
carried out by dividing the total evolution operator into short segments in such a way that 
the Hamiltonian does not change significantly. For instance, 

N-l 

Uit,to) = n (2-10) 

n=0 

where At = t/N. The spatial discretisation also has a direct influence on the time- 
propagation as the kinetic energy spectrum is different for different choices of Ar (see 
below). The most commonly used short- and long- time propagators are discussed briefly 
below, 

2.3.1 The second-order differencing scheme(SOD) 

In the second-order differencing(SOD) scheme[85, 88] which is also referred to as the ex- 
plicit method [188], the temporal evolution of the wavefunction is carried out as follows: 

^n+l ^ ^ ^(^3^ (2.11) 

Ti 

where and are the wavefunctions at (n-fl), (n— 1) and steps respectively. 

Therefore, once we know the wavefunction at (n — 1) and steps, we can use the above 
relation to evaluate the wavefunction at the (n-f-1)*^ step. Eq.(2.11) can be easily obtained 
by the following route: 






( 2 . 12 ) 
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where U and are the time-evolution operators. By expanding them using the Taylor 
series and truncating after the second-order term we get 

^n+l ^ ^ ( 2 . 13 ) 

a 

In order to use the SOD scheme, we need to know the wavefimction at (n — 1)^'* and 
steps [85, 88]. This scheme conserves both norm and energy. However, it requires a 
small time increment (At) for numerical stability as otherwise contributions from higher 
order terms in At will ficcumulate and make the scheme unstable. The stability criteria 
for step-size (At) can be obtained by considering the maximum energy {Emax) on the grid 
that has been introduced by the spatial discretisation. Emax is obtained by 


E, 


T = 

•^max 


= Tmax + Ki 


max max • ^ max 

,2 1.2 t 2 


(2.14) 


2m 


2m 


In the Fourier method 


TT 

Ar 


Substituting eq.(2.15) in (2.14) results in 


E = 

■^max — 




d-K. 


(2.15) 


(2.16) 


2m(Ar)2 

From time-energy uncertainty relation we can estimate the critical value of the time step 

h 


Atarit ^ 


E„ 


(2.17) 


Therefore, we have to choose At which is smaller than Afcrit- Any time step larger than 
this will make the scheme unstable. For all practical calculations Kosloff and Kosloff[85] 
recommended that 

Atcrif 


At 


(2.18) 
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be used. This method is easy to implement and it also conserves the norm of ^ and energy. 
It can also be applied to situations where the Hamiltonian is explicitly time-dependent. 
Another important point to be stressed is that the Hamiltonian must be strictly Hermitian. 
Therefore, it cannot be used in conjunction with absorbing boundary conditions (see below). 
In such situations one can use first order propagation scheme[91]. 

2.3.2 The split operator method(SOP) 

Feit and Fleck [82] introduced this scheme to solve the paraxial wave equation in opti- 
cal fibres. Subsequently they extended this scheme to solve the Schrodinger equation as 
well[83, 84]. According to this scheme, the time-evolution operator is approximated as, 

exp(—iAtH/ti) = exp{—iAtT/2h)exp{—iAtV/h)exp{—iAtT/2h) + O(At^) (2.19) 


Therefore, a single step propagation in time consists of the following steps: (i) the coordi- 
nate space wavefunction is Fourier transformed into momentum space and then multiplied 
at each grid point by half kinetic energy propagator, ea!:p(^^^), and then the wave- 
function is inverse Fourier transformed into coordinate space, (ii) potential energy part is 
then evaluated by multiplying ^ at each grid point by exp(^^^^) and (iii) step (i) is again 
carried out to complete a full step time-evolution. 

As a result of the noncommutability of the kinetic and potential energy operators, the 
error accumulates in third order of At. This propagation scheme is unitary and conserves 
only the norm. The critical time-step needed for time-evolution can be determined by 


Atcrit < 


TT 

3Ay 


where AV' = Vmax — Vmin- Vmax and Vmin are the maximum and minimum potential energy 
represented on the grid respectively. In comparison with the SOD scheme, we can choose 
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slightly larger time step in the SOP scheme. Also this method(SOP) is stable when an 
absorbing boundary is included. 

Since the Fourier transformation is carried out four times for each step propagation, 
this scheme is computationally more time consuming, especially for larger dimensions. 
However, in practice, except for the first and the last step, the half kinetic energy operator 
is merged so that the action of the kinetic energy operator is evaluated in full time step. 

2.3.3 The Chebyschev method 


The propagation schemes considered thus far are applicable only in situations where the 
dynamics takes place in short time scale. However, the Chebyschev propagation scheme 
is a global one, meaning thereby that the time step can be taken larger[91, 186, 189]. 
Therefore, the Chebyschev method is, in principle, useful for investigating long time 
dynamical events such as reactive scattering involving complex formation, multiphoton 
absorption/dissociation, vibrational predissociation etc. According to this method, the 
time-evolution operator is approximated in terms of a Chebyschev polynomial expansion. 
Since the Chebyschev polynomial is bound in the interval (-1,1), the Hamiltonian has to 
be renormalised by shifting its eigenvalues to this range. This can be obtained from a 
rough estimate of the eigenvalues of the Hamiltonian. For instance, the maximum energy 
represented on the grid in the FFT method is 


E„ 


2/i(Ar)2 


+ Vm 


( 2 . 20 ) 


and the minimum energy is 


Emin ^min 


( 2 . 21 ) 



16 


Thus, the average energy is 


E 


Emax + Ermn 


( 2 . 22 ) 


and the range is 


AE 


Efnax Exnin 


(2.23) 


The renonnahsation of the Hamiltonian operator is then obtained by 

H-IE 


Hr. 


AE 


(2.24) 


where I is the identity operator. Now the time-evolution operator can be rewritten as 


exp 


■iHE 


exp 


-itE^ 


exp 


'-iAEtHn 


(2.25) 


Expanding the second exponential term in the above equation in terms of complex 
Chebyschev polynomials leads to 


exp 


-iHr 


exp p Cr^Jr. r„ 


(2.26) 


where C„ = 1 for n = 0 and C„ = 2 for n > 1. The J„’s are Bessel functions of the first 
kind, of order n. r„(— i.H„orm) are complex polynomials satisfying the recurrence relation 


<l>n+\ = —2iHncrrm(f>n + <f>n+l (2.27) 

where = T„(— ijHn(?rm)^(0) with the condition <f>i = ^(0) and <l >2 = — *-ffnorm^(0). 
The Bessel function Jn{Q), where a = decays exponentially for n > a . Especially 
for larger a, Jn(cx) decays rather rapidly. Therefore, one has to include a few extra terms 
above n = a for convergence. 

For all practical purposes the minimum niunber of expansion terms should be about 40, 
as otherwise additional terms above n = a needed for convergence becomes dominant and 
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the method become inefficient [87, 186]. In the Chebyschev scheme, there is no restriction 
on the size of the time step. One can even complete the entire time-propagation in a single 
step. However, a major disadvantage with this global propagation scheme is that the 
intermediate results are not available. For example, processes such as photodissociation 
and resonance Raman scattering etc., mainly relay on the short time dynamics of the 
wavefunction on the excited electronic state. In order to apply the Chebyschev method 
to such short time dynamical processes we have to consider the time dependence of the 
expansion coefficients in eq.(2.26) as they are time-dependent. 

Using the time evolution operator obtained in eq.(2.26) , we can write the time evolving 
wavefunction as 

^{t) = (2.28) 

For processes such as photoabsorption, the cross section is computed by Fourier transform- 
ing the autocorrelation function (see below) as 

(t{u}) cx f e^('J'(0)|//l'f(t))dt (2.29) 

00 


Substituting for ^(t) from eq.(2.28) and carrying out the Fourier transformation leads to 

'h{uj — Wo) 


<t(w) oc ^ fe„ 

n=0 


AE 




(2.30) 


where wq = {E+ Eghift) /^- Eghift is the energy difference between the groimd and the 
excited electronic state and = $„(— iH„orm)’J'(0). fe„’s are the Fourier transform of the 
Bessel functions (J„). 

For n = even. 


— Wo) 

AE 




Wo) A 


- 1/2 


AE ) 


cos 


n sin 


-1 


^h{u) — Wq) 

AE 
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for n = odd, 


K 


'h{uj — ujo) 

AE 


= 4 1- 


'h{u — Wo) 
, AE~ 


- 1/2 


sin n sin 


-1 


'h{u — Wo) 

AE 


and 6o = 2 ^1 — ^ • Simulations carried out for the impulsive absorption 

spectrum of Csl yielded converged results [51]. Similar study has also been carried out for 
resonance Raman scattering as well [87]. This method is particularly efficient for calculating 
the absorption spectra of more complicated system. 

A major inadequacy of this approach is that it cannot be used for explicitly time- 
dependent Hamiltonians, as the evolution operator then cannot be expanded in Taylor 
series. 


2.3.4 The short iterative Lanczos(SIL) scheme 

This method consists of two steps: (i) the Hamiltonian operator is reduced to a tridiagonal 
matrix form and (ii) the numerical solution of the TDSE is obtained by diagonalising 
the tridiagonal matrix[91, 190-192]. The tridiagonal matrix is much smaller compared 
to the actual Hamiltonian matrix and it spans only a reduced subspace of the original 
Hamiltonian. The tridiagonal matrix is constructed via the following steps: 

H(f>o = oco<j)o -f- (2.31) 


and 


+ aj<l>j 4- ,j = 1, Nl-i (2.32) 

where Nl is the number of Lanczos recursions and (j){0) = \^(0). The are constructed 
such that they are orthonormal and the coefficients aj and /3j are obtained by the usual 
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Hilbert space inner products: 

aj = 

The reduced Hamiltonian operator has the following tridiagonal matrix form: 




f ao /3o 0 ••• 

A A 0 

0 /3i 02 A 

0 0 0 ••• 

VO 0 0 ••• 


0 \ 
0 
0 

OiNi-2 0 Nl-2 

0Nl-2 J 


(2.33) 

(2.34) 


(2.35) 


which is Nl x Ni compared to N y. N (where N is the total number of grid points) of the 
original Hamiltonian. 

In terms of we can write the evolution operator as 


U{At) = 


(2.36) 


For At of the order of tenth of a fe, 5-10 iterations are sufficient to construct the tridiagonal 
matrix. Diagonalisation of the evolution operator can be efficiently carried out so that the 
evolution operator becomes 

U{At) = (2.37) 

where Z is a unitary matrix that diagonalises and is the Hermitian transpose of Z. 

is the diagonal matrix of eigenvalues. Thus the time-evolving wavefunction is written 
as 

Nl-1 

^(t + At) = Z^e-^^^^^L/^Z(j}n 

n=0 


(2.38) 
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The number of recursions depends on the time step as well as the initial choice for <f>o - a 
good initial guess for leads to very fast convergence scheme. SIL scheme is uncondition- 
ally stable and it preserves norm and energy. Unlike Chebyschev method where the error 
is distributed uniformly, in this method the error is larger for larger eigenvalues. If the 
nmnber of recursions increases then there can be rapid loss of orthogonality of the recur- 
sion vectors. Therefore, it is always desirable to keep the recursion number small. Also the 
error starts building up, if the recursion vectors span outside the tridiagonal space. The 
error analysis carried out by Park and Light [94, 193] predicts that the Lanczos method can 
be used efficiently provided that a good initial guess is used to initiate the recursion and 
small enough time steps are employed so that the convergence becomes fast and storage 
requirement becomes small. 

2.3.5 The integral equation approach 

Kouri and coworkers [194-197] have pioneered the integral equation approach for solving 
the TDSE. An attractive feature of this method is that it is applicable for both time- 
dependent and time-independent Hamiltonians. In the following, we siunmarise briefly 
the methodology and the computational algorit hm of the integral equation approach. 

Let us consider a situation where the potential energy part is explicitly time-dependent 
as it happens in the interaction of radiation with atoms and molecules (see below). Then 
the TDSE will have the form 

(2.39) 

where “IXE.) is , the kinetic energy operator. For mathematical simplicity let us rewrite 
eq.(2.39) as 


(2.40) 
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By left-multiplying both the sides by we get 

giTtW I *) ^ Lf(R)yS{R, f) I = (2.41) 

However, we know that 

^ (i!, ()) = j ®(ii. t) + 

= I — ^ ■ ’ - + t)| (2.42) 

Substituting eq.(2.42) in eq.(2.41) results in 

(2.43) 

Upon integrating this equation between the limits t and f" leads to 

(2.44) 

Now, operating on both sides by yields 


f (2.45) 

/i Ji 

This equation is the central equation for the time evolution in the integral equation ap>- 
proach. In order to make this equation computationally feasible, we make use of the 
trapezoidal rule for the integral that is appearing in eq.(2.45). Thus, 






t 

h 


(2.46) 
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By substituting t" = t + t, - we obtain 
^{R,t + T) = 

- ^ i) + t + T)^(i?, * + r)) (2.47) 

On rearranging eq.(2.47) we get 

or 

<S(R,t + T) = (l + |^VtK,t + T)) e-<«''W'‘(l-^{t-R.«))®W«) (2.48) 

Thus the kinetic energy operator is flanked by the time- dependent potential energy op- 
erators. This approach is also referred to as kinetic referenced modified Cayley (KRMC) 
approach. It is important to note that the time-ordering property has been retained in the 
eq.(2.48). Since the terms -f + rfj and ^1 — ^i\R,t)j are the first two terms 

of the Taylor series expansion of exp{iTV(R,t+T) /2h) and exp{—iTV[R,t)/2h) respectively, 
we can replace them by their appropriate exponentials. Hence, eq.(2.48) becomes 

— I'T ^ “'2T' — I'T 

'^(R^t + r) = exp{—V{R,t + T))exp{—^R))exp{—V{R,t + r)MR,t)i2A9) 

This equation is similar to that discussed in the split-operator method except for the time- 
dependent potential energy part. The time-evolution can be accomplished by evaluating 
the kinetic and potential energy operators in the representations in which they are local. 

In a recent study, Ttuong et al. [92] have compared the efficiency and accuracy of the 
second and third order split-operator method (see below), the integral equation approach 
(KRMC) and the Chebyschev polynomial expansion method. Their study shows that the 
KRMC method is highly recommendable for short time dynamics as (i) this approach 
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requires only 2 FFT calls for each time step propagation, (ii) the trapezoidal rule that 
is employed in approximating the integral is a linear function of time; hence, one can 
choose larger time step as well, (iii) the Hamiltonian may be either time-dependent or 
time-independent and (iv) the norm and energy are conserved. However, for long time 
dynamics the Chebyschev method is preferable in comparison with the KRMC method. 
But if the potential under consideration is time-dependent, then one has to use the integral 
equation approach. 


2.3.6 The exponential split operator method 

Bandrauk and coworkers [181, 198, 199] have developed the exponential split operator 
method, as an extension of the split operator method discussed earlier(section 2.3.2). It 
is superior to its predecessor in that it is applicable for both time-dependent and time- 
independent Hamiltonians. 

For time-dependent Hamiltonians, the time-evolution operator will have the form 


, ^ /_j ft+At ^ \ 

U(t + At,t) = rexp(-^jf H{t)dt\ 


\njt ■ ■ / 

where Tis the time-ordering operator. By expanding eq.(2.50) one obtains 


(2.50) 


U{t -h At, t) = exp H{t)dt\ + Dz{t){iAt/Kf + Di{t){iAt/Kf -f . . (2.51) 

where Dz{t) and Di{t) are higher order terms in commutators of H{t). The unitarity 
property of the time-evolution operator allows us to omit the even order error terms in 
eq.(2.51), as they are real quantities and vanish. Therefore, lJ{t -t- At,t) becomes 

U{t + At, t) = ( -^ H{t)dt] + Dzit){iAt/hf + D 5 {t){iAt/hY + . . (2.52) 

In order to develop the new propagation scheme, let us split the time interval (t -f- At) 
into {t+sAt,t),{t+{l — s)At,t + sAt) and {t+At,t+{l — s)At) subintervals. U{t + At,t) 
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will then become 

Uit + At, t) = U{t + sAt, t)U{t + (1 - s)At, t + sAt)U{t + At,t+{1- s) AtX2.53) 


The time evolution operator for each of the subintervals can be obtained by making 
use of eq.(2.52) as 

U{t + sAt, t) = exp f — jf H{t)dt j + Dsit) x {sAtf + 0(Af)® 

= S2{t sAt,t) Dsit)] X (sAt)^ + 0(Ai)® (2.54) 


where Sz{t + sAt,t) = exp(—'isAtA/2h)exp{—iBlh)exp{—isAtA/2h) with A = Tl(i2) and 
Similarly, 


U(t "f" (1 — s') At, t “h sAt') = Sz(t “I" (1 — s)At, t -|- At) 

+[C' 3 (f) + A(*)] X [(1 - 2s)sAtf + 0(Af)® (2.55) 


where 


S 3 {t + (1 — s)At, t + sAt) = exp(— i(l — 2s)AtA/2h)exp{—iB/h) 

exp{—i{l — 2s)AtA/2h) 

with A = ilR) and B = t)dt and 

U{t + At,t + {1 - s)At) = 53(t + At,t+(l-s)At) 

+[C3(t) + i?3(t)] X (sAt)^ + 0{At)^ (2.56) 

where 


S 3 {t + At,t+ (1 — s)At) = exp{—isAtA/2h)exp{—iB/h)exp{—isAtA/2h) 
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with A = 3t-R) -S = Therefore, combining eq.(2.54), (2.55) and 

(2.56) we obtain 

i7(A<) = S 3 (sAt)S 3 ((l - 2s)At)S3isAt) + [25^ + (1 - 2sf] x [C 3 {t) + D 3 {t)]At^ 

+0(At)® (2.57) 

where C 3 {t) and D 3 {t) are higher order commutators. When s satisfies the equation 2s^ + 
(1 — 2s)^ = 0, the leading error in Ai^ vanishes even for a time-dependent potential. Thus 
the eq.(2.55) can be used to accomplish the time propagation. An important point to note 
is that the time-dependent potential which is appearing in equations (2.54-2.56) can be 
solved by using the middle point formula[199]: 

jf f{t)dt = Atf{t + At/2) + 0{Af) (2.58) 

More recently, Bandrauk and Shen[181] have applied this method to investigate the pho- 
todissociation of H 2 molecule by an intense laser pulse. The kinetic energy operator has 
been evaluated using the FFT method. Since the higher order split operator is unitary, 
the norm of the wavefimction is conserved for all the time steps. An added advantage of 
the higher order splitting scheme is that the phase of the wavefunction is also reproduced. 
Therefore, this is an ideal scheme for evaluating the eigenvalues. Also this method can be 
applied to studying nonadiabatic processes. 

In addition to these methods, recently, Fujimura and coworkers [200] have introduced a 
sliced propagator method for solving the TDSE with explicitly time-dependent potentials. 
According to this propagation scheme, the time evolution operator is sliced into several 
small segments in such a way that the time evolving wave packet becomes 


<f>{t + €) = U{€)(j){t) 

= exp{—iH{t)e/K)(j){t) 


(2.59) 
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where c is an infinitesimally small time step. The sliced-propagator, U{e), is expanded in 
Taylor series and the terms up to third-order are retained. Therefore, 


trie) = 1-F 


-iHit)€\ 1 (-iHit)e\ 1 (-iH{t)e 


h 


+ 


2! 


+ 


3! 


(2.60) 


This method was employed to investigate the above-threshold dissociation of Hj molecular 
ions. It was used in conjunction with FFT method for spatial evolution and the norm was 
found to be conserved for substantially small time steps. 


2.3.7 Long time dynamics 

In the earlier sections we have discussed various short time and long time propagation 
methods. However, if the process of interest is one of the following kind; (i) reactions 
involving intermediate complex formation (ii) nonadiabatic coupling between electronic 
states occurring at large internuclear distances (iii) unimolecular dissociation processes 
such as vibrational predissociation and (iv) systems having heavier masses, then one is 
expected to time evolve the wavefunction long enough to complete the process. One 
of the major difficulties with the long time dynamical processes is that a portion of the 
wavefunction (fast moving component) reaches the grid edge (asymptotic limit) while other 
portions of it are residing in the interaction region[201, 202]. Those components which are 
reaching the grid edge will show spurious oscillations and lead to error in the calculation. 
One way of circumventing this problem is to use larger spatial grid. However, this is not 
advisable as this makes the process computer intensive. Several authors have proposed 
various numerical methods (described below) to study the long time dynamical processes 
without extending the grid size. 

For instance, Neuhauser and Baer[203] have successfully demonstrated in collinear H 
+ H 2 reaction that a negative imaginary potential(NIP) can be used to absorb the product 
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channel wave packets whenever they are reaching near the grid edge. They have also shown 
that the inclusion of a NIP does not alter the reaction probability results. Subsequently, 
Neuhauser et al. [204-206] have explicitly derived the functional form of the NIP and the 
selection criterion for it. The widely used absorbing boundary has the form: 


/ Robs < R< Robs + ^Rabs 

I 0 R ^ Robs 


(2.61) 


where Vio is the height of the absorbing potential and AR is its width. Rj is the internuclear 
distance at which the absorbing potential starts building up. Invariably Rj is selected far 
away from the interaction region. 

The height and width of the absorbing potential are selected based on the energy(i^ 
of the time-evolving wave packet. It has been shown by Neuhauser et al.[204, 205] that 
the height of the absorbing potential must satisfy the following condition: 

(2.62) 


where fi is the reduced mass of the system. Similarly the width of the absorbing potential 
must be such that 

< AR 

where Emax ” 2/x(AR)^ ^max and Ejnin ~ 

Since the product channel wave packet is absorbed near the grid boundary, one has 
to use alternative methods to compute the final state distribution. Recently, Child has 
carried out a detailed analysis of the properties of the complex absorbing potential [207]. 
Several reactive scattering[94, 208, 209] and photodissociation[125, 170] calculations have 
been performed using the absorbing potentials and they have shown that the absorbing 
potentials are highly useful when the dynamics requires long time propagation. More 





28 


/ 

recently, Vibok and Balint-Kiurti have discussed various functional forms of the absorbing 
potentials and their applications[210, 211]. 

One of the major drawbacks in using the absorbing potential is that the information 
contained in the wavefunction is lost once it is absorbed. In order to circumvent this 
problem Heather and Metiu had proposed the wavefunction splitting algorithm[212, 213). 
According to this algorithm, the time evolving wavefunction is split into an interaction 
part, <), and an asymptotic part, before it reaches the grid boundary. The 

splitting is carried out efficiently using a splitting function (/(i?)) as follows: 

= fiRrr,)^{Rm,ti),m = l,2,...N (2.64) 

and 


= (2.65) 


where the subscript 1 refers to the splitting performed at time fi. The function f{R) is 
chosen such that it remains 1 throughout the interaction region and decays smoothly to 
zero near the asymptotic region. For instance. 


m) = 


( 2 . 66 ) 


1 + exp{—a(R — R)) 

The parameters a and R are chosen depending upon the size of the asymptotic region. It 
is worthwhile to add that the time evolving wavefunction is not altered by this splitting 
procediue as we have only multiplied ^{R,t) by 1 = f{R) + 1 — f{R). Since the time 
evolution operator is linear, both ^j{R,t) and ’J'^(i?,t) can be propagated independently. 
As we know that 'J'^(i?,f) is in the force-free region, i.e., V(R) = 0, its temporal evolution 
is obtained by multiplying it with the free particle propagator at each momentmn grid 
point as 

'$^^(#1 + nAt,pj) = exp '^AituPj) i = 1, 2, . . . JV (2.67) 



29 


where p/s are the momentum values discretised as pj = {j — N/2 — l)Ap,j = 1,2,...N 
with Ap = where L is the length of the grid in coordinate space. 

The time evolution of the interaction region wavefunction is carried out by applying any 
one of the time propagation methods discussed earlier. Further wavefunction splitting is 
carried out when a part of the time evolving wavefunction reaches near the grid boundary 
as 

^l,2{Rm,h) = = 1,2, ...TV (2.68) 

and 

'^A,2{Rm,h) = (1 — = 1, 2, . ..iV (2.69) 

• 

where <2 is the time at which the 2"^ splitting is performed. Again the asymptotic and 
interaction region wavefunctions are propagated in the same way as discussed above and 
the splitting is carried out whenever the interaction region wavefunction reaches the edge 
of the grid. The splitting frequency depends upon the size of the asymptotic region. 
If the asymptotic region is small, then greater number of splittings are required. This 
procedure is continued till all the wavefunction residing in the interaction region reaches 
the asymptotic region. In the end, after M splittings, we can reassemble the wavefunction 
in the momentum space in the following manner: 

nt/,Pi) = E j=h2,...N. (2.70) 

where the asymptotic pieces which have been separated at times fa, a = 1, 2, . . . M, are 
propagated from ta to tf by the application of a free particle propagator in the momentum 
space, and are added coherently. This method is easy to implement and it does not require 
much nximerical effort. 
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Metiu and coworkers have extended this method to determine the final state distri- 
bution in the case of polyatomic photo£ragmentation[212-214]. For further applications 
please see reference [94]. 

Thus far we have concentrated mainly on the elimination of the wavefunction reflec- 
tion fi'om the grid edge. The spreading of the wavefunction all over the grid is an inherent 
property of the Schrodinger picture wavefunction as it contains both kinematical and dy- 
namical information. However, if the TDQM calculations are performed in the interaction 
picture formalism, the wavefunction spreading can be avoided to a large extent. Since the 
kinematic information is removed in the interaction picture wavefunction (see below), it 
does not spread or translate in the coordinate space. As a result, one can chose smaller 
spatial grid and larger time step in the interaction picture formalism. Recently, Zhang 
and coworkers have utilised this method to study the vibrational predissociation[215— 218] 
and molecule surface scattering processes [219]. The interaction picture wavefunction is 
obtained from the Schrodinger picture wavefunction through unitary transformation as 

^i(t) = exp (iHot/fij (2-71) 

= exp exp '^s{0) (2.72) 

where H = Hq + V and Hq and V are the free part of the Hamiltonian and the interaction 
potential respectively. From eq.(2.71) it is clear that ^j{t) is obtained by back propagating 
'3?s(t) with the free Hamiltonian Hq. Interestingly, both the Schrodinger picture and the 
interaction picture coincide at time t = 0. The Hamiltonian in the interaction picture is 
obtained by 

Hi = exp(iHot/h)Vexp[-iHot/h) (2.73) 

It is important to note that Hj is explicitly time-dependent, even when V is time- 
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independent. If we differentiate eq.(2.71) with respect to t and multiply by ih, we obtain 
a Schrodinger type equation 

== mi{t) (2.74) 

This equation can also be written as 

= exp[iHot/h)i^st (2.75) 

Since the interaction picture equations (2.74) and (2.75) are similar to the TDSE eq.(2.1), 
the time propagation methods discussed earlier can be easily implemented. However, care 
must be taken, since the interaction picture Hamiltonian is explicitly time-dependent. 

Zhang and coworkers[219, 220] and Tannor and coworkers [22 1-223] have extended this 
method to obtain product vib-rotational state distributions in vibrational predissociation 
and photo dissociation processes. More recently, Williams et al.[222, 223] have developed a 
modified version of the above method known as nested-interaction representation and suc- 
cessfully apphed it for studying three dimensional (3D) (J=0) photodissociation dynamics 
of triatomic molecules. 

Apart from these numerical methods, more recently, Guo[59] has introduced a novel 
method while investigating the 3D (J=0) photodissociation of CH3L According to this 
method, the translational coordinate(i?) is shifted during the time evolution depending 
upon the location of the center of the wave packet, (R). Since the wave packet remains 
compact in most of the direct photodissociation processes, the center of the wave padcet 
can be estimated periodically. Whenever the center of the wave packet is larger than the 
center of the R grid, the grid will be shifted by an amount AR so as to acconnnodate the 
change. By shifting the grid in this way, the momentmn space picture is not affected as 
the total number of grid points are not changed. The norm has been monitored carefully 
and the grid edge reflections are successfully avoided by this technique. 



32 


2.4 Interaction of radiation with matter 

Interaction of radiation with matter is one of the fundamental problems in physics. It has 
been studied extensively by both semiclassical and quantum mechanical approaches. In 
the former, the atom or the molecule (matter) is treated quantum mechanically and the 
electromagnetic radiation is treated classically. Such a treatment gives accurate results 
for processes like photoabsorption and emission. But it is not adequate to describe purely 
quantum mechanical phenomena such as spontaneous emission or quantum beats etc. In 
the following sections we will show how the semiclassical and the quantum mechanical 
treatments yield the same formula for the single photon absorption process. Such a treat- 
ment is essential to see the correspondence between classical and quantum approaches and 
to obtain an insight into the dynamics. 

2.4.1 Semiclassical treatment 

As we stated above, the system (i.e., atom or molecule) is described quantum mechanically 
and the radiation is described classically. To start with, let us consider the relativistically 
invariant form of the Hamiltonian for the motion of a charged particle in an electromagnetic 
field[224, 225]: 

-ihV--A f -h e(i> (2.76) 

2m c 

where m is the mass of the particle with charge e and c is the velocity of light. A and <f) 
are the vector and the scalar potentials associated with the field. 

The Hamiltonian in eq.(2.76) can be rewritten as. 
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Prom the commutator relation of [A, -i/iV] we can show that, 


iW.(AV>) = ih{V.A)^ + A.iihVxp) 


(2.78) 


However, in Coulomb gauge we have V.A = 0 and <f> = 0. Therefore, 

ifiV.A = ihA.V (2.79) 


Substituting this result in eq.(2,77) gives. 


H = 


ihe . 




• (2.80) 


The dynamics of any system in the presence of the electromagnetic field can be obtained 
by solving the time-dependent Schrodinger equation (TDSE): 

■ft-Lp = + ^A.V+ + nr)} I ®W ) (2.81) 

where | ’3P(#) ) represents the state of the system at any time t. Here V(r) is the system 
potential. 

The vector potential A describes the time-dependence of the radiation field and can 
be defined as [224], 

A(r, t) = -h (2.82) 


where Aq is related to both intensity and polarisation of the radiation. The electric and 
magnetic fields of the radiation are related to A as follows; 


e(r,t) = i6i;{Aoe^’^*’ — A^e^ Tk.r+»art)j^ 

B(r,f) = (i/AXo){Aoe(*-^-’‘^) - (2.83) 


B(r,t) is always oriented perpendicular to €(r,t). 
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For a weak electromagnetic field we can neglect the term and therefore eq.(2.81) 
becomes 

The general solution of eq.(2.84) can be written as, 

I ^(^, ^) ) = E c-tW I Mr) ) (2.85) 

k 

where summation is over all the discrete and continuous sets of | ) and Ek refers to the 

energy eigenvalue of the field free Hamiltonian Hq for the system. 

For an electronic transition from a discrete level | V’i) to a final state | V’/)j the wave- 
function for the system is written as a linear combination of the initial and final states: 

I ^(r, t) ) = CiWe-®''* I V-iM ) + C/We-®/-/" I i>,(r) ) (2.86) 

Substitution of the above expression for |^(r,t)) in eq.(2.84) and further simplification 
yields, 

I Ur) ) + | V-/M ) 

= H(t)ci(t)e-®*/‘ I Ur) ) + \ ^,(r) ) (2.87) 

Now, left-multiplying with {tpf{r) | on both sides of eq.(2.87) and integrating over r, 
we obtain 

( 2 . 88 ) 

= ( V-y I H{t) I V’i ) ft(t)e-®'/* + { V, I H{t) I V/) c,(t)e-®'"’' 

The first term in the left hand side of the above equation vanishes, because | V’i ) and ] ^}>f ) 
are orthonormal. 
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For a system in the state |^i) before the field is turned "on", Ci{t = — oo) = 1 and 
Cf(t = — oo) = 0. In addition, for a weak field Ci(f) ~ 1 and Cf(t) « 1. Under these 
conditions, eq.(2.88) becomes[3]. 


dcf{t) 

dt 

Substituting for ^(t) we obtain. 




(2.89) 


dcf(t) 

dt 


n me 

4(^/1 V’i ) 


me 


me 

e 

me 


+~ 


me 


{ xPf I V | V-i ) 

{ I Aoe'^^'^ V i 

( ijf I A^e-*^-'’. V I tpi ) e^(Ej+^-Ei))t/h 




(2.90) 


where 

0‘li = — {V'/|Aoe*'.V|*i) 

me 

0 ‘fi = — {V>/|Ao*e-‘KV|Vi) 

■' me 

Therefore, 


Cf{t) = afij" + 

_ 2 'j ( J 1 

-i{uj - w/i) j “H i(w + lOfi) J 


(2.91) 


(2.92) 


where u/i = (Ef — Ei)/h 

The first term in eq.(2.92) vanishes when Ef = Ei + hu, i.e., for absorption process and 
the second term vanishes when Ef = Ei — hu, i.e., for emission process. Therefore, the 
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first and the second tema in eq.(2.92) represent the absorption and the emission processes 
respectively. 

Photoabsorption 

The probability for photoabsorption from state |^i} to state | V’/) is, therefore, given 
by 


Pfiii) 






{U - UfiY 


2 fl — + 1 


, ,, f 2 — 2cos(lu) — w 

M 1 > • sa " 


(W - Ufi}^ 


,2 [ 1 - 

= 1— 


- Ufi]t) 


Making use of the fact [226] that the 


lim — = 7d6[uj - Ufi) 

t-*oo — Ufiyt 


eq.(2.93) becomes. 


Pfi{t) = 2Trt\afi\H{u-Ufi) 


Therefore the transtion rate 


Rli = 


(2.93) 


(2.94) 


(2.95) 


(2.96) 


The exponential part in the expression for o/j in eq.(2.91) can be expanded in Taylor 
series: 


„ak.r _ ^ j. ^ 


(2.97) 
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Optical transitions typically involve atomic waveftmctions which extend over the order 
of the first Bohr radius of the atom, i.e., about 1 A. Since the wavelength (A) associated 
with optical transitions are of the order of 1000 A, k (= 27r/A) will be of the order of 10® 
cm"^ Thus the quantity k.r will be very small for r < 1 A. Therefore, we can replace e**' *’ 
by unity, i.e., all other terms including the first order term in eq.(2.40) are neglected. This 
approximation is known as, "electric-dipole approximation" or "long wavelength approxi- 
mation". But this approximation becomes less accurate as the frequency of the radiation 
increases [224, 227]. Therefore, within the electric-dipole approximation O/j becomes, 

Oifi = ^{%pf\Ao.V\rpi) (2.98) 

Replacing V by ^f, (viz., p = m f = — i hV), the above equation becomes, 

TUC fl 

= ^(^/l Ao-f IV’i) (2-99) 

However, in the Heisenberg equation of motion[224] both f and r are related through, 

r = ^[r.Ho! (2.100) 

Thus, 

1 ^ 

a/i = 

ocfi = ^w/i(V’/|Ao.r|V’t) (2.101) 

Substituting for a/i in eq.(2.96) we obtain 


( 2 . 102 ) 
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The energy removed from the radiation field of frequency, w, will be [228], 


^rad 


huRfi 

j wjj ( xjjf I Ao.r I ipi ) \^6{ijj - uji) 


(2.103) 


Now, the photoabsorption cross section, cr(a;), can be defined as the rate of removal of 
energy from the radiation field due to transition of interest divided by the intensity of the 
incident radiation [38], 


cr(w) 


Erad 

7(w) 

;^ 27 r j wJilAopI {‘ipf I e.r | - Ufi) 


(2.104) 


where Ao = € is tiie polarization vector of the radiation. However, from the Poynting 

vector theorem [229] we know 


|Ao|^ = (2.105) 

Substitution for |AoP in eq.(2.105) gives, 

«^(‘^) = ^^^l(^/|e.(er)lV’i)P^(w-a;/i) (2.106) 

We can also replace the term, er, by the electric-dipole operator(^) and ^ by the vacutun 
permittivity constant(eo). Thus, the most general form of cr(a;) is written as, 


( il’i I c-A* I V’i ) - t^/i) (2.107) 

2.4.2 Quantum mechanical treatment 

The total Hamiltonian of the system in the presence of radiation is given by[227, 230] 

H = Hrr^ + Hr + Hi 


(2.108) 
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where is the Hamiltonian of the atom/ molecule in the absence of any external field, 
Hr is the Hamiltonian of the radiation field, 

Hr = (2.109) 

with a^ and a the creation and the annihilation operator respectively. Hi is the Hamiltonian 
for the radiation-matter interaction and within the electric-dipole approximation Hi is 
given by, 

/ * \ 1/2 

m) = ( 2 , 110 ) 

where V is the volume of the cavity in which the radiation field is considered and eo is the 
vacuum permittivity constant. Photoabsorption causes transition from state \'4>i] Hi) to 
lipf^No), where | tpi ) and Itpf) are the initial and the final state of the molecule and the 
states I Ni ) and | No ) are the photon eigenstates before and after absorption. 

The time-development of the system in the presence of radiation can be obtained by 
solving the TDSE, 

= {So + «(*)} !*(*)> (2.111) 

where Ho = Hm + Hr. Since Ho and Hi{t) do not commute with each other the usually 
stated solution 

= U{t, 0)\^i0)) = 1 'i'(O) > (2.112) 

becomes inapplicable in the present case. However, we can use the time-ordered integral 
equation form of the time-evolution operator [231, 232]. Thus, 

U{t,0) = Uo{t,0)-^ f dt'Uo{t,t')H{t')Uoit',0) 

n Jo 

t dt' t dt"tro(t,t')^(*')t^o(t',t")^(^")^o(t", 0) -t- . . . (2.113) 

h Jo Jo 
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Alternatively, 

U{t,0) = ^ (2.114) 

n Jo 

t f + • • • 

fx Jo Jo 

where the subscript i referes to the matter-field interaction Hamiltonian. For transitions 
between two orthonormal states |^i) and IV*/)? ^be zero order term in eq.(2.115) does 
not contribute towards the transition probability. However, the higher order perturbation 
terms do contribute significantly. For a weak radiation field we need to concentrate only 
on the first-order perturbation term. Therefore, the transition amplitude is given by, 

= (2.115) 

h Jo 

Substituting for in the above gives, 

= ( V-/; I - ^ 1^1 ' ii.e{a - a*) j | jv, ) 

= ( r dt' { V>/; No I e-‘^»‘/V.e(a - a^) | V’ii N > e^(E,+Eo-[Ei+E,+nu]y /n 
\ V J Jo 

= ( V-/; No I e-‘^»‘/V.e(a - a^) | n ) jT (2.116) 

where Ei and Ef are the initial and the final eigenvalue corresponding to states | -tpi ) and 
|'0/) and jSo + aw and Eq are the eigenvalues corresponding to the photon states \Ni) 
and I No ). The latter is considered as a vacuum state. Fmther simplification of eq.(2.117) 
leads to, 

( . j \ 1/2 * f p-i(oJ-u;fi)t 1 'J 

— ) ( No I e-«</V,(a - a*) I 1(2.117) 
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Upon integrating over radiation field variables we obtain, 

The probability of transition at any time t is, 


(2.118) 


Pd*) = |W<j>(‘)P 




(w - UfiY 


In the limit < — > oo. 


(2.119) 


HmP/i(t) = (^■^~^\{‘tpf\fi.e\'ipi)f27rt6{u;-u}fi) (2.120) 


In the limit f — > oo, the last term in the right hand side of eq. (2.120) reduces to a 6-function 
representation as in the semiclassical description. The rate of transition will be given by, 



dPfi 

dt 




( 2 . 121 ) 


Now we can define the photodissociation cross section, u'(w), as the rate of removal of 
energy fi:om the radiation field of frequency a;, divided by the intensity of the incident 
radiation[38]. 


(t(w) 


hujRfi 

Nihujc/V 




Nitujc/V 

^ K V'/ 1 M-e \^i)f - t^/i) 


( 2 . 122 ) 


(2.123) 


The transition matrix element fifi = {ipf\ti.e\ipi) plays an important role in deciding 
whether the given transition is allowed or not. fifi also depends on the angle between the 
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vectors // and e, i.e., 

/t.e = |/i||€|cos0 

|(V’/|/i.e|V'i)r = l(V'/lMe|V’i)l^cos^^ (2.124) 

If the radiation under consideration is not polarised, then e will be oriented randomly. 
Therefore, cos^6 can be replaced by its average, 1/3 [224]. The photoabsorption cross 
section(for the unpolarised light) can thus be rewritten as, 

= 3^^l(V’/l/ie|V'i>p6(a;-w/i) (2.125) 

This expression is widely used to compute the photoabsorption cross section over all 
photon frequencies (cj) [3, 225]. The cr(cj) obtained by eq.(2.126) is the partial photoabsorp- 
tion cross section for transition into state /. The total photoabsorption cross section is 
calculated by summing over all the final states, 

<^i^) = w/i) ( 2 . 126 ) 

Nearly two decades ago, Heller and coworkers[75] had shown that a(a;) could be com- 
puted easily in the time-domain. Substituting the integral form of the ^-function, 

6{u}-Ufi) = r (2.127) 

2i7r J—oo 

in eq. (2.127) yields: 


(2.128) 
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where E = Ei and | ) = e I )• Hamiltonian corresponding to the 

final state /. co has been replaced by l/47r. \(j>i) = /xe | ^*1 ) 1® frequently referred to as the 
"promoted state". The final state wavefunction is normalised such that Y,f\i^f) (V’/l = 1- 
It can be seen from eq.(2.129) that the photodissociation cross section is the "full" Fourier 
transform of the autocorrelation function (C(t) = {<j>i\<j>i{'t)))- Therefore, from a practical 
point of view, <t{uj), can be computed readily by computing C{t) and obtaining its Fourier 
transform. 

2.4.3 Nature of the wavefunction prepared on the excited state 
under the influence of light 

In this section we will consider an electronic transition form the ground electronic state to 
an excited electronic state. Within the electric-dipole approximation the transition-dipole 
moment is given by /i,e(f), where e(f) is the electric field vector defined as. 


e(t) = A^ef{t)cos{u£) 


where Aq = 2 Aq. f{t) defines the shape and duration of the laser pulse. The TDSE of the 
radiation-coupled system is given by [233], 


ih 


dt 


V’cW 


At.e(f) He ) V Mi) ) 


(2.129) 


where the subscripts g and e refer to the ground and the excited electronic state respec- 
tively. We can rewrite the TDSE as, 


where, ^(f) 


= {Ho + H{t)}ni) 

eft 



( H, 0 

V 0 He 


and H{t) = 


0 

/i.e(f) 0 


(2.130) 
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Eq.(2.131) is similar to eq.(2.112), for which the solution has been obtciined already. 
Since it is fair to assume that the laser pulse begins to interact with the molecule only at 
time t = 0, there is no wavefunction on the excited state at t = 0, i.e., ‘ipe{0) = 0. Therefore, 
the time-evolving wavefunctions on the ground and the excited electronic state[131, 234] 
respectively are: 

f *'£-“•<‘-•'>'"^((>-'^•*'^>5(0) (2.131) 

rt */0 

We have assumed a weak field and hence the higher order terms have been neglected in 
writing the above equation. On substituting for in eq.(2.132) we get, 

n Jo 
/i */0 

= -^Ao (2.132) 

n Jo 

where E = Eg + huj. In order to arrive at eq. (2.133) we have made use of the rotating wave 
approximation (RWA)[227], cos(wt) = This approximation is valid near resonance 

condition {hu = Ef — Ei). From eq.(2.133) we can draw the following conclusions: the 
ground state wavefunction, “ipgiO), is constantly transferred to the excited state as long 
as the light is ’on’. While the wavefunction is being transferred, the initially arrived 
wavefunctions are propagated under the influence of the excited state Hamiltonian. 

Since in practice excitation is caused by pulsed lasers or continuous wave (CW) lasers, 
we consider two limits: (i) 5-function and (ii) continuous wave (CW). 

5-functioii limit 

Experimentally it is not possible to obtain a 5-function laser pulse. But theoretically it 
is important to consider this limit as it helps us understanding the Heller’s approach[75], 
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which has been employed extensively in the present work. In the 6-function limit, when 
t'=0 the laser pulse suddenly promotes the ground state wavefunction to the excited state 
and hence eq.(2.133) will be modified as, 

= (2.133) 

where ^^(0) = fi. e^,(0). 

CW limit 

In this limit, the Ught is on for ever, i.e., t — > oo and Aq = f{t) = 1. Therefore, we 
have, 

V»e(t) = lim (2.134) 

t—>oo fl Jq 

By substituting u = t — t' in eq.(2.135) we obtain, 

= lim r due’^“/''e-‘^‘“/Va(0) 

t^ooh Jo ‘ 

= lun^e-’^/^jT (2.135) 

Here <f>g{u) = e“’^‘“/^^j(0), the time-evolving wavefunction on the excited electronic 

state. Eq.(2.136) suggests that the ground state wavefunction is suddenly transferred 
to the excited electronic state as in the 6-function case, it is then propagated on the 
excited state and after substantial time-evolution, the stored wavefunction is half Fovuier 
transformed into energy domain. This energy domain wavefunction was previously defined 
by Heller and coworkers as the Raman wavefunction [72, 234-236], 

ME) = jT 


(2.136) 
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Thus we can rewrite the wavefunction prepeired on the excited state by a CW laser as, 

V’eC*) = (2.137) 

t — >00 

We can see from the above equation that the Raman wavefunction is what is prepared 
by the CW laser on the excited state. Now let us turn our attention to the photoabsorption 
cross section defined in eq.(2.129), 

On further simplification we obtain, 

= (2.139) 

Eq. (2.140) clearly demonstrates that the Heller’s approach to determine the photoabsorp- 
tion cross section falls in the CW limit and hence it is applicable to most of the experimental 
situations. 


2.5 Dynamics involving more than one excited (elec- 
tronic) state 

2.5.1 Adiabatic and diabatic representations 

We have discussed above the computation of ( t { u }) for a direct photoabsorption/dissociation 
process. However, in several photodissociation processes one frequently encounters more 
than one electronic state participating in the dynamics[99, 103, 107]. In order to describe 
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the coupling between several electronic states, we have to go beyond the Born-Oppenheimer 
approximation. This can be achieved by going over to either the adiabatic (non-crossing) 
representation or the diabatic (crossing) representation. In the former the nuclear kinetic 
energy operator, Tn{R), is non-diagonal and the potential energy operator, V{R), is diago- 
nal. In contrast, in the diabatic representation Tn{R) is diagonal and l^/J) is non-diagonal. 
For a two state system, we can write the TDSE in the adiabatic representation as, 

.. a / Mt) '^r+Tu Tn 

) \2ml^ r,, + 

where Tu,Tu,T 2 i and T 22 are defined as follows: 

Ti, = 

Til = + 

and 

^22 = 


^l(^) 

1p2it) 


(2.140) 


Here -01(2) and 0“(2) nuclear and the adiabatic electronic eigenfunctions for the state 

1(2) respectively. and V 2 are the adiabatic potential-energy curves for states 1 and 2 
in that order. 

The perturbations arising from the off-diagonal matrix elements are referred to as non- 
adiabatic perturbations. It is worth emphasising that the off-diagonal matrix elements in 
this case are exclusively due to the nuclear kinetic energy operator. 

For the same two state system we can write the TDSE in the diabatic representation 


as, 

1 ( Mt) \_S( TiiR) 0 ^12 

aiV02(t)7~U 0 niR)) \ %2 


V’i(i) 
Mi) ) 


(2.141) 
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where Ti{R) and T 2 {R) are the kinetic energy operators for nuclear motion on surfaces 
1 and 2 respectively. Vn( 22 ) and ki 2 ( 2 i) are the matrix elements ^ | | 

( '^ 1 ( 2 ) I ^d I ^ 2 ( 1 ) )• ^ 1 ( 2 ) ^1(2) are the diabatic electronic and the nuclear wavefunctions. 

In this representation the off-diagonal matrix elements are due to the nuclear potential 
energy tenn(i.e., Hei)- Perturbations arising from these off-diagonal matrix elements are 
termed as electrostatic perturbations. In the adiabatic representation the potential-energy 
curves may approach each other closely but they will not cross (avoided crossing). But in 
the diabatic representation the PE curves do cross whenever degeneracies occur. 

The two representations are related to each other via the following unitary transfor- 
mation: 


4>i \ _ f cos9{R) —sin9{R) ^ 

^2 J \ sin9{R) cos9{R) y \ ^^2 


where the i?-dependent angle 9 is obtained from 


tan29{R) 


2 Vi 


12 


Pn-P22 

The adiabatic PE curves are related to the diabatic PE curves through 


ya ^ (Vii + Va)±V(bi-V22)H4V;| 

1,2 2 

2.5.2 Time-evolution on coupled electronic states 


(2.142) 


(2.143) 


(2.144) 


The various methods described in section 2.3 for time-evolving the wavefunction on a single 
PES can be extended readily to systems involving more than one PES which are coupled 
with each other. For a 2-level system, for instance, at any point on the spatial grid, the 
SOD scheme will be as follows: 

+ Vn ] ( Mt.) \ ,2 

V ; V tfe(<n-l) ) h \ fi(ie) + {'22 A J 
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where 1 and 2 refer to the two states and the rest of the variables have their usual meaning 
as in eq.(2.142). 

Different authors have used the different schemes: SOD, Lanczos, Chebyschev, etc., 
depending upon the nature of the non-adiabatic dynamics. In addition, Coalson and 
Kinsey[78] and Alvarellos and Metiu[123] have proposed alternative strategies to study 
non-adiabatic interactions. 

Coalson and Kinsey[78] proposed the wave packet perturbation theory for a two-state 
system and it can be extended to more complex systems. 

They partitioned the diabatic Hamiltonian that appeared in eq.(2.142) into the unper- 
turbed and the perturbation parts as, 


where Ho = 


Hn ^0 

0 H22 


H 

andH = I 


Ho + & (2.146) 

Here Hn and H 22 are the zeroth order 


0 Vu 

V 21 0 


Hamiltonians for states 1 and 2 respectively. Letting ^ ^ substituting 

it in eq.(2.142) yields 




(2.147) 


By making use of the integral equation form of the time-evolution operator (cf. eq.(2.115)) 
we obtain. 


h Jo 

~iw) lo ' lo ' 

^-m(t '-t /g-iFot "Aijr(o) + . . . ( 2 . 148 ) 


where the exponential operator is. 


/ Q-iHnt/h Q \ 

0 ) 


(2.149) 
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By substituting the respective terms in eq.(2.149) and assuming that at t=0 there was no 
population in state 2 (i.e., V’ 2 ( 0 ) = 0) we obtain: 

( 0 ) 

+ (2.150) 

and 

V'2(*) = + (2.151) 

Coalson and coworkers [121, 237] have successfully carried out several numerical studies 
by making use of Gaussian wave packet method. These studies showed that the WP 
perturbation theory approach was adequate only when the coupling between the potential 
energy curves was weak. It failed to describe strong coupling situations. Also the GWP 
remained Gaussian only on harmonic or quadratic potentials. One can apply this method 
to anharmonic potentials, only if the dynamics takes place in short-time scales. Therefore, 
if the non-adiabatic transitions occur slowly, then we cannot use the GWP method. A 
more detailed account of the strength and weakness of the GWP method for studying the 
Qon-adiabatic transitions is given by Ramakrishna[237]. 

Alvarellos and Metiu[123] suggested an algorithm which is an extension of the split 
operator method (see section 2.3.2) to study curve crossing problems. In their scheme, 
:;he exponential operator is calculated in the representation in which it is diagonal. For a 
;wo-state system, the algorithm consists of the following two steps. 

Step 1: The diabatic nuclear wavefunctions are first transformed into adiabatic wavefunc- 
:ions and then multiplied by the phase factor 


-iAtVif/a 


^-iAtV^lh 


Then the resultant 



wavefunctions are transformed back into the diabatic ones: 
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V XiiR) ) \ Q A j ^ '4,^{R) ) 


(2.152) 


where T and T ^ are the transformation matrices to carry out the transformation between 
diabatic-to-adiabatic and adiabatic-to-diabatic representations respectively. They are de- 
fined as 

rpi _ ( cosB{R) sin6{R) \ rp _ ( cos9{R) sin9{R) \ 

-sin9{R) cos9{R) j ” V -sine{R) cos91r) ) 


and x’s in eq.(2.153) are the diabatic nuclear wavefunctions before and after evaluating 
the potential energy part respectively. The i?-dependent angle 9 has been defined earlier. 
Step 2: The wavefunctions obtained at the end of step 1 were Fourier-transformed into the 
momentum domain and then the kinetic energy operation was carried out by multiplying 
at each grid points with ( ^ g-iAtRA^/2m j > where and k 22 are the wave 

numbers corresponding to states 1 and 2 respectively. An inverse Fourier transformation 
is then carried out to complete one step time-evolution, 


\^P2{R,At) ) \ 0 [x2iR) ) 


(2.153) 


The efficiency of this method will improve if the following conditions are met: 

(i) total number of grid points must be less and the grid increment (Ai?) may be large; 

(ii) the dynamics must take place in short-time scales. 

This method can easily be extended to more complex systems. The relevant adiabatic- 
to-diabatic transformation matrix elements are obtained from the diabatic potentials (Cf. 
equations 2.143-2.145). 

The major time consuming step in this approach is the transformation between the 
two representations. However, a new scheme proposed recently by Broeckhove et al[238] 
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makes use of Lie algebraic manipulation to determine the time-evolution without switching 
between the adiabatic and the diabatic representations. According to this approach the 
diabatic potential energy matrix given in eq.(2.142) can be written in terms of Pauli 
matrices as, 


+ V22)(To + Vi2<7i -f- 2 (^11 ~ Vl2)(^3 


(2.154) 


where the Pauli matrices are given by, 


0-0 


1 0 
0 1 


; = 


0 1 
1 0 


; <^2 


0 —i 

1 0 


; 


1 0 
0 -1 


For mathematical convenience we can use the following simplifications to write the diabatic 
potentials; 


1=0 

== -j- l'20'2 + I^3^3 


(2.155) 


where 


^0 = 

— Vl2 
1/2 = 0 

t^3 = \{Vn-V^) 


and (To -3 are the Pauli matrices. Once we have represented the diabatic potentials using 
the summation as shown in eq.(2.156), we can always use the following relation to compute 


the exponential operator, 

exp{i ^ Uiai) = exp{iuo) icrocos\u\ -h i (^UiO-i 

Vi=i ' 


»=o 


sin\v\ 

"bT 


(2.156) 


\v\ 


3 

I 

<t=l 




(2.157) 
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This equation allows to calculate the potential as well as the kinetic energy operator in 
the diabatic representation and more efficiently computes the wave packet time-evolution. 
Using model studies Broeckhove and coworkers [238] have shown several application of this 
method. This scheme has one major drawback that is this method is applicable only for 
a two state problem. For multielectronic state cases determining the exponential operator 
through eq.(2.156) will be difficult. 

2.6 Predissociation 

If the molecule is excited to a state which is quasi-bound with respect to dissociation, then 
it may undergo predissociation. In diatomic molecules the following two types of predis- 
sociation have been observed. They are (i) rotational predissociation and (ii) electronic 
predissociation[16, 116]. For a rotating molecule the inclusion of centrifugal potential re- 
sults in a centrifugal barrier to dissociation. Therefore the molecules are quasi-bound with 
respect to dissociation; they can undergo predissociation by tunneling through the barrier. 
However, in the case of electronic predissociation, a molecule which is residing in a bound 
state that is well below the dissociation limit can dissociate by mixing with a continuum 
level of another electronic state. The two electronic states involved in this predissociation 
may or may not cross with each other. 

Due to predissociation the spectral characteristics will be altered significantly. For 
example, the intensity of the spectral lines will be lowered, spectral lines may be diffused 
or the energy of the spectral lines may be shifted. Experimentally, in a collision-free 
environment, it is possible to dete rmin e the predissociation lifetime (7^) or the line-width 
(r, more precisely the full width at half-maximum (FWHM)). Most of the calculations 
till this date have been carried out using TIQM approach[116, 239, 240], the primary 
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reason being the is of the order of s. Naturally, the TDQM approach 

without approximation would require propagation for a long time. However, with the 
developments of newer algorithms (see section 2.3) and computer architecture, it is now 
possible to investigate such long-time dynamics using TDQM approach[139, 167]. 

Recently, Villarreal et al.[168] have proposed a wave packet golden rule treatment to 
determine the line- width of the vibrational predissociation of Ne...ICl. According to this 
approach, the time-independent Fermi-Wintzel golden rule expression is transformed into 
time-domain as has been done for the photoabsorption process (cf. section 2.4.2). The 
Fermi-Wintzel golden rule expression for the predissociation rate (k^) is. 


'f' = ^ m- Ei) (2.158) 


where \'tpi(E } ) and ^ are the initial bound and the final continuum state of tbe 

molecule respectively. Similarly Ei and Ef are the eigenvalues corresponding to the initial 
and the final state of the molecule respectively. Vij is the interaction potential correspond- 
ing to the coupling between the bound and the continuum PE curves. On replacing the 6 
function in eq.(2.159) by its integral form and simplifying it further we obtain. 


1 

kr = - (2.159) 

ft J —oo 


where (j>i{0) = and (j){{t) = Hf is the Hamiltonian for the final 

state. 

Similarly the predissociation line-width, T{E), can be determined by 


r{B) = 

where r(jE) and are related through 


1 


(2.160) 


r{E) = 


27rc7*^ 


(cm ^) 


(2.161) 
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In summary, computation of r(£') involves the following steps: 

1. the initial state eigenfunction (^j)is placed on the final continuum state after 
operating it with the coupling potential(t^j). 

2. the autocorrelation function, C(t) = (^i(0)|<;6f(t)), is computed at each time-step 
during the time-evolution. 

3. r(jE) is then obtained by Fourier transforming the C(t). Since we have chosen a 
specific initial state 'ipi, we have to extract r(jE/) corresponding to the eigenvalue 
Ei. 

Further details regarding the computation of T{E) are given in chapter 5 for hydroxyl 


radical predissociation. 



Chapter 3 


Nal Predissociation: A Case Study Using 
Femtosecond Laser Pulses 

3.1 Introduction 

Technological developments in recent years have enabled the study of chemical reactions in 
real-time using femtosecond laser pulses[29, 241]. For instance, Zewail and coworkers[242] 
have investigated several unimolecular and bimolecular reactions using ultrashort laser 
pulses and they have demonstrated that the transition from reagents to products could 
be monitored while the process is occurring. Among various femtosecond transition- 
state spectroscopy (FTS) studies, Nal predissociation is particularly well studied[243]. 
The importance of Nal predissociation emerges mainly due to the nature of its potential 
energy (PE) curves shown in figure 3.1 [244]. The ionic ground electronic state(X^Eo) 
and the covalent excited electronic state (A O''') exhibit a diabatic crossing near 7 A. In 
the adiabatic picture these two PE curves show an avoided crossing and the excited state 
becomes a (quasi)bound state. Therefore, if the laser pulse prepares a wave packet on the 
excited adiabatic state, it will undergo(damped) periodic oscillations. Depending upon 
the Landau-Zener curve crossing probability (Pi^), whenever the wave packet representing 
the state of the system reaches near the crossing point, some part of it leaks into the lower 
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adiabatic state and results in the neutral products Na and I. The remaining part of the 
wave packet will continue to oscillate in the adiabatic well with a characteristic oscillation 
period. Thas, ttie (electronic) predissociation of Nal has become an ideal candidate to 
investigate using the ultrashort (ptunp - probe) laser technique[242]. 

In the experimental studies, Zewail and coworkers employed a femtosecond pump laser 
pulse to prepare 'Nal on the covalent excited state and the system was later probed using 
another femtosecond laser pulse. When the wavelength of the latter was in resonance with 
the Na D-line, i.e. , A = 589 nm, free Na produced due to predissociation was monitored (on- 
resonance probing) , However, when the probe wavelength was either red or blue- shifted 
from 589run, the Na atoms which were trapped in the adiabatic well (i.e., the transition- 
state, [Na...I]^’*) were probed (off-resonance probing). The laser-induced fluorescence(LIF) 
signal corresponding to the off-resonance probing showed an oscillatory pattern persisting 
up to 10 ps and the on-resonance probing signal showed a stair-case like increment [243]. 
The oscillation period was found to be dependent on the photoexcitation wavelength. 
Subsequent experiments carried out by Zewail and coworkers [245] over photoexcitation 
wavelengths (pump pulse) in the range 300 - 390nm enabled them to characterise the 
excited state PE cmve and also to observe the dynamical evolution of the system in real- 
time. In addition, Cong et al.[246] showed that the damped oscillations of the off-resonance 
LIF signals were recurring every 30 ps (revival time) and were persisting for about 40 ps 
or so. 

Subsequent to these experimental investigations several theoretical studies have been 
carried out using classical [247] and quantum mechanical [138, 139, 248] approaches. Choi 
and Light[248] and Engel and Metiu[139] have carried out time-dependent quantum me- 
chanical calculations of Nal predissociation. The former authors assumed that the pump 



60 


wave packet (WP) prepared under various pump pulse excitation wavelengths. We have 
chosen the following photoexcitation wavelengths: 300, 321, 344, 360 and 390 nm. Since 
resonance excitation occurs at 320 nm, these wavelengths would cause Tiear, vevy Titav 
and far from resonance excitation. 

In the next section we describe the computational method employed in our study. 
Results obtained are discussed in section 3.3 and the conclusions are given in section 3.4. 

3.2 Computational method 

Our computational approach is similar to that of Metiu and coworkers [133, 139], The 
wave packet prepared on the excited electronic state by a pump laser pulse is given by (see 
section 2.4.3) 

v>2(*) = -V r (3.1) 

ft J—oo 

where subscripts 1 and 2 refer to the ionic ground state and the covalent excited state 
respectively. ipiiO) is the ground electronic state eigenfunction. We have taken a harmonic 
oscillator wavefunction corresponding to v" = 0. £'„"=o represents the zero-point-energy of 
the Nal molecule. We have not considered the vibrational excitation of the molecule in its 
ground electronic state as our interest is mainly in studying the time-evolving character- 
istics of the wave packet prepared by the laser pulse, fj, is the electric dipole vector of the 
molecule. e{t') is the electric field vector of the electromagnetic field represented by 

e(t') = (3.2) 

where sq relates to the intensity of the field and f{i/) defines the shape of the pulse, w/ is 
the firequency of the laser pulse. We have used the Gaussian pulse 


f(t') = 


(3.3) 
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where bj is related to the width of the pulse and to is the center of the pulse. We have 
used the same pulse parameters as used earlier by Engel and Metiu[139]. They are: bf 
= 1.1 X 10 ^ fs ^ and to = 80 fs. The operator is the time-evolution operator 

corresponding to the covalent excited state. 

As we have discussed earlier, the ionic ground state and the covalent excited state 
are nonadiabatically coupled to each other. We use the diabatic representation[238] to 
describe the time evolution on the coupled electronic states as it is computationally conve- 
nient. Since Choi and Light [248] have pointed out that the diabatic coupling between the 
electronic states are not important during the pulse preparation, we ignored them during 
the pulse preparation step. The Gaussian pulse considered in the present calculation dies 
off in about 160 fs. Therefore, we have performed the pulse preparation calculation up 
to 160 fs. After the completion of the pump pulse preparation, we have carried out the 
time-evolution of the WP on the coupled electronic state, where we have explicitly in- 
cluded the nonadiabatic coupling. We have employed the second-order diflferencing (SOD) 
scheme[85, 88] for time evolution. Since Nal is a heavy mass system, i.e., its reduced 
mass, (iNai = 19.5 amu, we had to choose a small grid spacing (Ai2). Our numerical grid 
consisted ofi? = 3.0 (0.03) 33.69 au (1024 grid points). We have carried out time-evolution 
for 3 ps using tim e steps of 0.53875 fs. In order to avoid WP reflection from the grid edge 
we have used a damping function as was suggested by Guo [250]. 

3.3 Results and discussion 

As has been mentioned in the Introduction, we have considered a single photon absorption 
from the ground electronic (ionic) state to the covalent excited electronic state. Fimctional 
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forms of the relevant PE curves and nonadiabatic coupling potentials used are[251]: 

V,(R) = + + 


-^1 , (^Na + o:j 2ajvi,a7 ' 


47reo [ R 




(R+R^y^ 

Vi2{R) = 




Relevant parameters have been taken from the work of Wang et al.[251] and they are 
reproduced in table 3.1. The PE curves and the coupling potentials are illustrated in figure 
3.1. Since we have no information about the transition dipole moment, we have taken it 
to be imity. In figure 3.1 we have indicated the points a-e referring to the photoexcitation 
wavelengths: 300, 321 , 344, 360 and 390 run respectively. It is clear from this figure that 
the latter two wavelengths correspond to off-resonance excitation. 

The norm of the wave packet during the pulse preparation is given in figure 3.2 as a 
function of time for different pump pulse wavelengths. Since the laser pulse is essentially 
off in about 160 fs, we have plotted the norm only up to this time. The variation of 
the norm for different photoexcitation wavelengths show the characteristic behaviour of 
near, very near and off -resonance photoexcitation. This trend was predicted earlier 
by Williams and Imre [131] for a model system and it can be understood in terms of the 
transients prepared by the laser pulse. In order to explain this transient formation, let us 
rewrite the integral given in eq.(3.1) as a siunmation: 

N ( 

M*) = E ( 

n=0 \ 

where (j){0) = ^^i(0) and the following substitutions are made for time: At = dt', t = NAt 
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Table 3.1: Parameters for the PE curves[251]. 


Ai 

2758.0 eV 

Bi 

3.2919 eV^/* au 

C] 

579.270 eV au® 

Di 

1.072 

Pi 

0.6809 au 

O^Na 

9.2789 au^ 

Oil 

36.7515 au^ 

^1 

5.312 eV 

A 2 

3150.0 eV 

B 2 

5.1212 eVi /12 

C 2 

45540.13 eV au® 

P2 

0.80824 au 


3.18 eV 

Ro 

1.1527 au 

hx 

0.055 eV 

A 2 

0.1920 au-2 

Rx 

13.782 au 


X 
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and t' = nAt. 

Eq.(3.7) can be rewritten as 


N 


^2{t] = Y,^2{nAt) 


n=0 


(3.8) 


where 


-iAt 


il)2{nAt) = I —j— I (3.9) 


Alternatively we can write eq.(3.8) as 


V’2(<) = '02(O)+^2(A<)+^2(2At)+...+V'2(AAt) 


(3.10) 


Thus it is clear from eq.(3.10) that the pump pulse prepared WP, ‘tpii't), contains infor- 
mation of all ^ 2 (^Af). Since we do not know at what instant the pump pulse photon is 
absorbed by the molecule, we add all the transient WPs produced on the excited elec- 
tronic state[133]. Under resonance or near resonance photoexcitation, i.e., when Xex = 
321 nm and 344 nm, the WPs which are arriving on the excited electronic state add up 
constructively and lead to a steady increase in population. This maximum population 
remains stable for all the time even when the pulse is off. However, in the off-resonance 
cases shown in figures 3.2a, 3.2c and 3.2d, the norm increases when the pump pulse reaches 
its maximum amplitude and then decays with time. This is because the WPs which are 
arriving after the pulse maximum is reached (t > to)) interfere with the components which 
arrived earlier. They are invariably out of phase with each other, and therefore they inter- 
fere destructively and the norm decreases[66, 131]. Nevertheless, it is interesting to note 
that some part of the WP survives firom this interference and continues to remain on the 
excited electronic state. Therefore, even when the pulse is off there is a small amount 
of norm stabilised on the excited electronic state. However, the norm showed in 3.2e for 
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390 nm photoexcitation, rises to a Knaxiiniiin at t = tg and falls down subsequently. This 
phenomenon is purely due to destructive interference. It is also important to notice the 
magnitude of the norm for all photoexcitation wavelengths. Also the interferences between 
WP components are particularly iooportant as it will affect the population significantly. 
This point will be clarified further (see below). 

After the pulse excitation step is ower, we have plotted the probabiUty density of the 
WP as a function of internuclear (distance for different photoexcitation wavelengths in 
figure 3.3. The magnitude of the probability density varies drastically with the excitation 
wavelength considered in this study. Another interesting feature of figure 3.3 is that the 
width of the probability density profile decreases as the excitation energy is decreased. 
This is because the number of superposition states contained in the WP prepared on the 
excited electronic state decreases ■when the photoexcitation energy is decreased. Plots of 
|‘02(-R)P for 360 and 390 nm, excitation in figure 3.3d and 3.3e respectively require some 
additional comments. They show that the probability density amplitude is particularly 
low. In addition, for the excitation wavelength 390 mn, there is an additional peak near R 
= 11.5 au. This unusual characteristic of the WP prepared by the laser pulse is purely due 
to the transient eflfect. As we hawe mentioned earlier, the early arriving components of 
the WP start leaving the Franck-Coadon region while other components are still arriving. 
Thus there are constructive and iestntctive interferences. Any amoimt of WP that 
has escaped from the Pranck-Coadon region starts moving towards larger internuclear 
separation. This happens especially when the photoexcitation wavelength is away firom 
resonance[131]. As a result of this tnnsient effect we see a decrease in the probability 
density amplitude as well as tanxisuaal preparation of the WP. When these WPs are time- 
evolved on the coupled electronic states (i.e., diabatically coupled excited and ground 
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electronic states) they exhibit some interesting features which are described below. 

Engel and Metiu had already pointed out that the coherently prepared wave packet 
('02(^)) remains intact even after long time[139]. Since the oscillating WP gets refocussed 
at the ionic and the covalent turning points, the WP remains Gaussian for a substantially 
long time. We have also observed this characteristic in our studies. However, we have 
observed some noticeable differences when we move away from resonance excitation. In 
figures 3.4a and 3.4b, we have illustrated the probability density of the time-evolving WP 
as a function of R corresponding to the photoexcitation wavelengths 300 mn and 390 
nm respectively. The upper and lower panels of figures 3.4a and 3.4b axe the probability 
density profiles of the WP on the covalent and the ionic state respectively. As we have 
mentioned earlier, the pump pulse prepares the WP in about 0.16 ps. Therefore, at t = 
0.16ps there is only the piunp pulse prepared WP on the covalent electronic state. As the 
time progress, the WP spreads and a part of it is transferred to the ionic state at aroimd 
0.26 ps. During the crossing, a small amount of WP is leaking and travelling towards larger 
R on the covalent excited state. Since its magnitude is relatively small, it is not visible 
in these plots. At a larger time, i.e., around 0.56 ps, the WP on the ionic state reaches 
the turning point. The WP becomes narrow at this time and then starts moving towards 
lower R. It reaches the crossing point around 1.16 ps. a part of the WP has already been 
transferred to the covalent state. After the WP reaches the covalent state, it moves towards 
lower R and collides with the repulsive wall and produces severe oscillations at 1.26 ps. 
Thus the intact WP fractiues into several small firagments. However, these fragments add 
up again and continues to oscillate between ionic and covalent turning points. Since the 
pump pulse prepares a coherent WP on the excited electronic state, it acts like a classical 
particle for a long time. 
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Figure 3.4a. Probability density of the time-evolving WP on both covalent and ionic 
states after the pulse preparation (photoexcitation wavelength corresponds to 300 nm). 
Upper and lower panels correspond to the WP on covalent and ionic states respectively. 
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Figure 3.4b. Probability density of the time-evolving WP on both covalent and ionic 
states after the pulse preparation (photo excitation wavelength corresponds to 390 nm). 
Upper and lower panels correspond to the WP on covalent and ionic states respectively. 
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This is, however, not valid when the photoexcitation wavelength is completely away 
from resonance ( particularly when detuning is done towards lower energy). For instance, 
the figures given in figures 3.4b for 390 nm show several interference patterns. In order 
to understand them let us focus our attention on the WP prepared by the pump pulse as 
shown in the first panel of figures 3.4b. It has two peaks and its origin has been explained 
earlier. The early arrived component of the WP ( marked as EA) has already reached near 
the crossing point while the late arrived one ( marked as LA) is in the Franck-Condon 
region. At 0.26 ps, the early arrived one is transferred to the ionic state, whereas the 
late arrived component is still remaining on the covalent state. Within the next 100 fe, 
i.e., 0.36 ps, some part of the late arrived WP is transferred to the ionic state and it 
moves towards the ionic turning point. However, the early arrived WP component had 
already reached the ionic turning point and started moving towards the covalent state. 
Around 0.46 ps, these two oppositely travelling WPs collide with each other and produce 
the interference. Due to this interference, the probability density amplitude decreases and 
the total WP fractures into several small WPs and time-evolves further on the ionic and 
covalent states. However, when they reach the repulsive wall of the covalent excited state, 
they do not add up constructively, instead they will further split into multiple WPs as 
shown in figure 3.4b for times 0.86-1.06 ps. For 360 run pulse excitation also we have 
observed the WP spreading in about 0.5 ps. Therefore, based on our studies we can 
draw the following conclusion, when the photoexcitation energy is much lower than the 
resonance excitation energy, the WP prepared by the laser pulse spreads in a short while 
unlike the one prepared near the resonance excitation energy. Both the transient effect 
and the pump pulse energy are responsible for the observed behaviour of the WP. 

We have carried out the time-evolution of the WP up to 3 ps. We have determined the 
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population of the Na which is trapped in the adiabatic well using the formula, 

rRx 

Pboujidit) — / lV’2(-R)t)pdi? (3-11) 

where Rrnin Rx aic the minimuni R value represented on the grid and the covalent- 
ionic crossing point respectively. Pi^nd{t) values obtained in this calculation are plotted as 
a function of time in figure 3.5a-e. For 300, 321 and 344 nm excitation( figures 3.5a, 3.5b 
and 3.5c), Pbmndi't) obtained in the present calculation shows oscillation periods which are 
in good agreement with the previous results[29, 139]. However, the bound Na populations 
shown in figures 3.5a— c, differ firom the experimental results. The experimental peaks 
are broader and the oscillation periods are slightly larger than ours. However, a better 
agreement can be obtained if we include the Boltzmann average of the populations. Since 
the WP prepared by low energy excitation spreads immediately during the time-evplution, 
we notice a difference in the Pboundi^) behaviour as a function of time for 360 and 390 nm. 
These plots shown in figures 3.5d and 3.5e are not in accord with experimental results 
of Zewail and coworkers [243]. In order to simulate their results one has to perform a 
two-photon calculation. 

3.4 Conclusion 

We have carried out TDQM calculations within the first-order perturbation theory and 
computed the pump pulse prepared WP on the excited electronic state corresponding to 
the photoexcitation wavelengths: 300, 321, 344, 360 and 390 nm. These photoexcitation 
wavelengths are near, very near and far away from resonance excitation wavelength (i.e,, 
320 nm). Therefore the WP prepared by laser pulses of different wavelength exhibits 
different characteristics. Norm of the WP prepared during the pulse preparation, time- 
evolution characteristics of the WP and bound Na population have been computed. Our 
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Figure 3.5. Bound Na population (Pi^cmnd) a function of time, (a), (b), (c), (d) and (e) 
correspond to photo excitation wavelengths 300, 321, 344, 360 and 390 nm respecti vely. 
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results are interpreted in 
obtained from the bound 
results. 


terms of the transient effect. In general, the oscillation period 
Na population for 300, 321 and 344 nm agree well with earlier 



Chapter 4 


Channel Control in a Chemical Reaction. Vibrational 
Enhancement of I*/I Branching Ratio in HI 

Photodissociation 

4.1 Introduction 

Several strategies have been proposed over the years to control the outcome of a collision 
event, to effect one channel preferentially over the other. Considering one of the simplest 


collisions, that between an atom(j4) and a diatom(RC) 

are possible depending upon energy and other factors: 

as an example, several outcomes 

A + BC-^ A + BC 

(NR) 

-^AB + C 

(Rl) 

-^AC+B 

(R2) 

— > A B C 

(cro) 


The simplest means of channel control known to chemists for a long time has been 
to increase the relative translational energy (Etrans) of reactants by supplying heat. For 
Etrans < Eth where Eth is the threshold energy for the exchange reaction, the only channel 
available is NR. For Etrans > Eth the exchange channel R1 opens up. Similar constraints 
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hold for R2 and CID channels. Depending upon the system, charge transfer and other 
channels open up at different Etrans- In the last 2-3 decades it has been realised that 
for many reactions the reagent vibrational excitation enhances the reaction cross section 
more selectively than E'trans [252]. The reagent rotation plays a subtler role[253]. For 
many reactions, rotational excitation of a reactant molecule leads initially to a decrease 
in the reaction cross section. Invariably this trend is followed by rotational enhancement, 
becoming comparable to vibrational enhancement. 

While such strategies rely upon initial state preparation of reactants, substantial effort 
has gone into examining the possibility of controlling the collisional outcome by intervening 
at the transition state[254]. We must add that the success via this route has been rather 
limited so far. 

In the case of unimolecular reactions, one is always concerned about the possibility of 
energy randomisation. Attempts have been made to understand the factors responsible 
for intramolecular vibrational relaxation (IVR) and to look for means of attaining mode 
selectivity [255]. 

Several schemes have been proposed in the last few years to achieve channel control. 
Rice and coworkers[117, 256] proposed a pump and dump scheme which envisages elec- 
tronically exciting a molecule (using a pump laser), letting the system time-evolve on the 
excited state and stimulating de-excitation (dump) onto the lower channel. Using a model 
system they showed (numerically) that channel control could be achieved by adjusting the 
time delay between pump nd dump pulses. Gadzuk[257] pointed out how the same could 
be accomplished in an adsorbate molecule by altering the workfuncion of a metal surface by 
depositing alkali metal atoms, without using laser pulses. Briuner and Shapiro[156] pointed 
out how it might be possible to utilise the coherence property of a laser for the same pur- 
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pose. Several workers ( for example, see ref. [155, 258, 259]) have examined the possibility 
of using differently shaped pulses to attain the same goal. More recently, Rabitz [260] and 
coworkers have examined the possibility of inverse quantum mechanical control. That is, 
depending upon the choice of products, work out (inversely) the requirement of the shape 
of the field and then put it to practice. 

Attempts have also been made to obtain one set of products preferentially over another 
by combining initial state preparation with photodissociation. For example, using model 
calculations, Imr e and coworkers[50, 52] pointed out the predominance of one isotopic 
channel over the other in the photodissociation of HOD, arising from initial vibrational 
excitation. Grim and coworkers [261] have demonstrated experimentally vibrationally me- 
diated photodissociation in HOD leading to a higher OD/OH branching ratio. Shapiro[99] 
had pointed out the importance of including non-adiabatic coupling while investigating 
the photodissociation of CH 3 I. Schatz and coworkers [53, 56, 262] showed that, vibrational 
excitation of the molecule in its groimd electronic state could influence the branching ra- 
tio, r, for F and I formation, in the photodissociation of CH 3 I. Kash and Butler [263] 
demonstrated that vibrational excitation prior to photodissociation serves as a probe of 
the nature of non-adiabatic coupling. 

Photodissociation of hydrogen iodide (HI) in the near UV region has been studied ex- 
tensively by theory[265-268] and experiment [269-273]. It is of interest in both dynamical 
and spectroscopic studies. In dynamical studies, for instance, the translationally hot hy- 
drogen atoms produced in the HI photodissociation are used as a source for the H atom to 
study H atom exchange reactions [273, 274] and translational -> vibrational energy transfer 
processes [275, 276]. 

As early as 1930, Mulliken[264, 265] predicted that the photodissociation of HI involved 
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photoabsorption from the ground electronic state (X ^Eq) to excited electronic states: 

and ^JIq. Among the excited states, photoexcitation to and states led to 
H( S) + I( P 3 / 2 ) to Hq state led to H(^S) + I*(^Pjy 2 )- He also clarified that the 
electronic transitions ^Eq ^IIi, ^Hj as perpendicular transitions and ^Eo — » ^Hq as 
a parallel transition. Therefore, according to this interpretation H + I* products are 
essentially produced from the parallel transition. Subsequently various experimental and 
theoretical studies have been carried out to understand the photodissociation mechanism 
and to obtain the nature of the potential energy curves for the system. 

For example, Ogilvie employed a 6 - function approximation to the spectral profile 
and constructed a single repulsive excited electronic state[266]. Clear et al.[270] applied 
photofragmentation spectroscopy and estimated that 36 ± 5 % of I were produced in the 
photodissociation of HI at the excitation wavelength, Aei = 266.2 nm. Also the mea- 
sured kinetic energy distribution and the product angular distribution were found to be in 
agreement with Mulliken’s interpretation. Schmiedl et al.[271] implemented laser Doppler 
spectroscopic technique and obtained similar results. The 1* /I branching ratio (T) ob- 
tained by them at Ae^ = 266 nm was ~ 1.5. The photofragmentation mechanism proposed 
in the former studies was in accord with the earlier predictions of Mulliken. However, 
the photofragmentation spectroscopic studies carried out by Van Veen et al.[272] at Xex 
= 193 nm, 222 nm and 248 nm, suggested a slightly different mechanism. Accordingly, 
the observed photodissociation of HI was ascribed to the following electronic transitions: 
^Eo ^Hi, ^Eo -> ^Hi, ^Eo ^Ho and ^Eo ^Ej. The former two transitions were 
of perpendicular type and led to H -f- I formation and the later two transitions were of 
perpendicular and parallel types respectively and they both led to H -f- 1* formation. All 
these excited states, except ^Hoj were nonadiabatically coupled to each other. Since the 
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coupling potentials were weak, the photodissociation cross sections and the I* /I branching 
ratio calculated by Van Veen et al. without including nonadiabatic couplings were foimd 
to be in good agreement with experimental results. Mulliken also considered the -» 
3 Ei transition in his studies[264]. However, he assumed that the state was lying higher 
in energy than the other states. Recently, Wittig and coworkers [100] have carried out high 
resolution spectroscopy studies on the photofragmentation of HI and HBr and the results 
obtained by them are in good agreement with the photofragmentation spectroscopy mea- 
surements of Van Veen et al. Therefore, it can be concluded the location of the ^Ei state 
and its nonadiabatic coupling with other electronic states are playing vital rule in deciding 
the outcome of the photodissociation. 

By making use of the experimental values such as the photodissociation cross sec- 
tions and I* /I branching ratios at different excitation energies. Levy and Shapiro[267] 
constructed semiempirical potential energy curves and nonadiabatic couplings for the dif- 
ferent states of HI. Subsequently they employed the TIQM approach and obtained the 
photodissociation cross sections, the T/I branching ratio and the anisotropy paxameter(/3) 
over the energy range 30,000-50,000 cm“^. The results obtained by them were in good 
agreement with the earlier studies. An important conclusion from their study was that the 
inclusion of nonadiabatic coupling was highly essential to estimate T accurately. Once T as 
a function of the excitation energy is known, one can select Aex and produce H atoms with 
specific translational energy for use in dynamical experiments. In addition, understanding 
the mechanism of branching might help us to look for some means to control T . 

Since the previous study carried out by Levy and Shapiro[267] was restricted only for 
the ground vibrational state (v” = 0) of HI in its ground electronic state, we decided to 
investigate the HI photodissociation in detail, including the effect of vibrational excitation 
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of HI on the nonadiabatic coupling using the TDQM approach. Since the method of 
computing the photodissociation cross section has been discussed explicitly in chapter 2 
(section 2.4.2.), we only state the salient features and present the results and discussion in 
the next section. 

4.2 Results and discussion 


As mentioned in the Introduction HI photodissociation involves photoexcitation from the 
ground electronic state (^Eo) to the excited states ^Hi, ^Hi, and ^Hq. The potential 
energy curves are shown in figure 4.1. Except for the ^Ho state all the other states are 
nonadiabatically coupled to each other. As we have discussed in the methodology chap- 
ter(see section 2.5), the multielectronic state dynamics can be easily solved in the diabatic 
representation. Therefore, we have carried out our calculations in this representation. 

The functional form of the diabatic potential energy curves for the excited states and 
the nonadiabatic coupling potentials axe taken fi:om Levy and Shapiro’s work. They are, 

mm) = AOexp{-B{0{R-R.)) + c{0 (4.1) 



where = ^ni,^ni,^Ei and ^Hq. The parameters used to fit the potential energy 
ctuwes are given in table 4.1. 

The ground electronic state of HI is represented by a Morse curve and the groxmd 
vibrational state (v”=0) wavefunction is transposed on the excited state potential energy 
curves given by Levy and Shapiro [267] and time-evolved using the corresponding excited 
state Hamiltonian. We have made use of the transition dipole moments reported by the 
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Table 4.1: HI potential energy parameters 
Ground electronic state parameters: 

De, = 0.117 450 259, 0g = 0.926 910 853, = 3.040 872 9 

Excited electronic state parameters: 



"Hi 

'Hi 



A(0 



0.04174 

0.04664 

m 



1.588 

1.843 

m 



0.034 650 9 

0.034 650 9 

MO 



0.009 983 

0.011 451 
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same authors. The nonadiabatic effects are fully taken into account by solving the multi- 
electronic state time-dependent Schrodinger equation: 


ih 


dt 


N 


^2 


W3 


^4 ) 
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Hn Vi2 Vi3 
^1 H22 V22 

V 0 0 


33 


0 
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\ 



to 



/ 



(4.3) 


where the subscripts 1,2,3 and 4 represent 3ni,ini,3Ei and states respectively. The 
second derivative of the wavefunction with respect to i?, the internuclear distance of HI, was 
obtained by the FFT method[85] and the time-evolution was accomplished by the second 


order finite difference(SOD) method[88]. The spatial grid consisted of = 0.1(0.1)25.6 
au and the time-evolution involved a total of 27,384 steps with each time step(Ai) = 


0.00218fs. 


The time-evolution was followed graphically for 11,000 steps by which time the pop- 
ulation in each state, reflected in the value of levelled off as shown in figure 4.2. 

This ensured that the effect of diabatic crossings on the wavefunction was fully taken into 
account as has been done by others[126]. For the wavefunction at the 11,000th time step, 
the autocorrelation function ('$'(T)|'$'(T’-|-t)) was evaluated for the remaining 16,384 (=2^^) 
steps. Here T refers to ll,000Ai. In order to compute the partial photodissociation cross 
section, the choice of T is irrelevant. But for computing the F, T must be greater than 
r, the time taken by the wavefunction to traverse the diabatic crossings [126]. The partial 
photodissociation cross section for each electronic channel was obtained by Fourier trans- 
forming the respective autocorrelation function and the results are shown labelled 1,2,3 
and 4 in figure 4.3a. All these results are in excellent agreement with those obtained by 
the time-independent route by Levy and Shapiro[267]. The value of F was obtained by 
taking the ratio of the sum of the partial photodissociation cross sections [79] for I* and 



lation(X1 0 






(a) v'=0 j ] |(c) v"=2 ] I (d) v^' = 3 | [(e) v"=4 



states respectively and the subscript v” indicates the vibrational state of HI. Results of 
branching ratio T = (r{H + r)/(T{H + i) for different v” are included in the top panels. 
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I channels and the results included in figure 4.3a are also in accord with the previously 
reported values. It is worthwhile to mention that due to the presence of nonadiabatic 
couplings the intensity of the IIi and are affected appreciably (intensity borrowing) 
though ’III is not affected much. 

The anisotropy parameter (^) is an important experimental quantity which charac- 
terises the nature of tfie electronic transitions[18, 128, 272, 277]. In the center of mass 
frame /3 parameter is related to the angular distribution of the products through 

[1 + /9P2 (cos^)] (4.4) 

where 0 is the angle between the detection direction and the direction of the polarisation 
vector of the incident light (e), (Tq is the total cross section, 47r is a normalisation factor and 
P 2 is the second order Legendre polynomial. In the case of a purely parallel transition /3 
= 1 and the net electronic angular momentum change during the electronic transition will 
be An = 0. On the other hand for a pmely perpendicular transition = -0.5 and the net 
electronic angular momentum change will be AO = ±1. Therefore, in HI photodissociation 
Ho transition is of purely parallel type whereas IIi,^ Hi and ^Ej transitions 

are of purely perpendicular type. For a parallel type transition the intensity is proportional 
to cos^Q. But for a perpendicular type transition it is proportional to sin^0. K a particular 
product channel originates from parallel and perpendicular transitions as it happens in the 
H -f- r channel, then we can write the j3 parameter as 

/3 =A - 0.5 B (4.5) 


where A and B are the percentages of parallel and perpendicular transitions respectively. 
Similarly, if /3 parameter is known A and B can be obtained from it by 


A = 


( 2 ^+ 1 ) 


B = 


(2 - 2/3) 


3 


5 


3 
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By making use of eq.(4.5), we computed the /3 parameter over the energy range 34,000 
- 52,000 cm-’. The results are shown in figure 4.4 along with the experimental values 
obtained by Van Veen et al.[272]. Our results are in good agreement with the experimental 
values[272] and the time-independent quantal results[267]. In the low energy region T 
formation is predominantly due to parallel transition whereas in the high energy region it 
is due to perpendicular transition. 

We have repeated the above calculations for different initial vibrational states of HI 
and the results are included in figure 4.3b-4.3e for v”=l-4. In general, one finds that 
there is an order of magnitude increase in T with increase in v”. Further, for each v” there 
are v”-|-l peaks in the T versus energy plot. For v”=l, the peak at the higher energy is 
particularly strong. For v”=2 and 3, this peak becomes less intense but the peaks at the 
lower energy become stronger. For v”=4, the higher energy peak is the strongest with 
two additional peaks being strong and two weak. Much of the structure in plots of r(F?) 
can be understood from the structures in plots of i=l, 2,3,4. For convenience, we have 
grouped the results for i=l and 2 corresponding to the H -|- 1 channel in one part of the 
figure and for i=3 and 4 corresponding to H -f I* in another. The major features in the cr*„ 
curves can be understood in terms of the "reflection principle" which predicts v"-t-l peaks. 
The additional structures in the cr’n curves can be traced to the nonadiabatic coupling 
between the different electronic states. 

In order to further elucidate the mechanism underlying the vibrational enhancement of 
r we have plotted in figure 4.5, the steady state population, as determined fi'om the 
value for distances beyond the diabatic crossings, on the four excited states as a function 
of v” . It is clear from this figure that the population on the state increases gradually 
with increase in vibrational excitation. However, on the ^Hi and ’Hj states the population 
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shows an oscillatory behaviour as a function of v«. These results in turn can be understood 
by considering the nature of the initial vibrational wavefunctions prepared on the excited 
state ^Ei, for example, for different v” as illustrated in figure 4.6. Particularly from the 
wavefunction depicted for v”=4 in figure 4.6e, it is clear that a significant portion of it is 
present beyond the crossing point. Since the ^Ei and states cross (diabatically) at R 
= 3.8 au, vibrational excitation of HI prior to photodissociation assists the wavefunction 
in escaping from the non-adiabatic coupling region. Therefore it is understandable that 
r(E) is strongly dependent on v”. 

4.3 Conclusion 

The time-dependent quantum mechanical calculation carried out for HI photodissociation 
from the ground vibrational (v"=0) state of the ground electronic state shows a good 
agreement with time-independent quantum mechanical results. The results obtained for 
the photodissociation of vibrationally excited HI shows an overall increase in P with v” 
and that it is highly structured as a function of the photoexcitation energy. These results 
have profound implications in that the branching ratio for two competing channels in a 
photodissociation process can be controlled by vibrationeil excitation. Also, this study 
suggests a means of enhancing the production of H atoms with a particular recoil energy 
for use in dynamics experiments. 
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Figure 4.6. Illustartion of the initial vibrational wavefnnctions transposed on the excited 
state vis-a-vis the non-adiabatic coupling region for different v”. For the sake of clarity, 
we have included only the potential-energy curves for the ground electronic (^Eo) &iid ITi 
and states. 







Chapter 5 


Photodissociation and Predissociation Processes in 

Hydroxyl Radical 


5.1 Introduction 

Hydroxyl radical(OH) has been observed frequently in cometary gases and in interstellar 
clouds. It is one of the major constituents in the earth’s atmosphere and in interstel- 
lar media. Also it is one of the important sources for 0 and H atoms in astrophysical 
environment[278, 279]. Therefore, understanding various dissociation pathways in OH 
becomes particularly important. Photodissociation and photon-induced predissociation 
are the major dissociative processes that occxir in terrestrial atmosphere. Studies along 
this direction were initiated more than forty years ago[280]. However, precise informa- 
tion regarding the potential-energy (PE) curves for the system has become available only 
recently. The important PE curves involved in the photodissociation and predissociation 
processes in OH are shown in figure 5.1. Ab initio calculations carried out by Dalgarno and 
coworkers [282-284] showed that photoabsorption from the groimd electronic state (X H) 
to the excited states, l^E”, l^A, and 2^11, led to direct photodissociation. Pho- 

todissociation cross sections(cr) computed by them were later shown to be a factor of two 
smaller than the experimental results[284, 285]. 


92 




R(au) 


Figure 5.1. Potential-energy curves involved in photodissociation and predissociation 
processes in OH. 
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Predissociation in OH takes place[278] in the first absorption band {A^I.+-X^U) due 
to spin-orbit coupling between the bound state and the repulsive states ^E", l^E" 
and ^n, shown in figure 5.1 as dashed lines. The (A-X) band predissociation has been 
studied quite exteasively by theory [287-289] as well as experiment [290-298]. For instance, 
Gaydon and Wolfhard[280] observed an anomalous decrease in spectral line intensities 
corresponding to all rotational levels of v’ = 2 and higher rotational levels of v’ = 1 while 
investigating the radiative lifetimes of the - X^II) transition and they attributed 

it to predissociation. Subsequent experimental and theoretical studies also confirmed the 
occurrence of a weak predissociation near v’ < 3. In addition to this weak predissociation, 
Carlone and Dalby[289] and Czarny and coworkers [291] observed a strong predissociation 
near v’ > 5. Czarny and coworkers suggested that the observed predissociation was due to 
spin-orbit interaction and they also made an estimate of the coupling between the A^E'*' 
and the ^IT state. They had also estimated the predissociation lifetimes to be of the order 
of 10'^' sec. Based on experimental observations [290-295] and ab initio calculations of 
Michels and Harris [286], it was proposed that the weak and the strong predissociations 
in OH were due to spin-orbit coupling between the bound state A^E'*' and the repulsive 
states ^E~ and ^H respectively. 

Sink et al.[287] employed a multichannel scattering calculation to determine the de- 
pendence of predissociation linewidths (F) on the rotational-vibrational levels of A^E"*". 
They included the ab initio values of the spin-orbit coupling between A^E'*' and ^E , 
^E“ and ^H. However, their spin-orbit coupling potentials were not explicitly dependent 
on the OH bond distance. Recently, Yarkony[298] has carried out an extensive configu- 
ration interaction (Cl) calculation to determine the spin-orbit coupUng matrix elements 
and their influence on the predissociation lifetimes. Unfortunately, the results reported 
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by him were restricted to the state and that too for the first four vibrational levels 
of the state. Since (A - X) band laser-induced-fluorescence is used extensively to 

monitor OH produced in collision-free environments [299], there has been a lot of interest 
in understanding the details of the predissociation. 

We decided to investigate both photodissociation and predissociation processes in OH 
using the TDQM approach to obtain an independent estimate of ( 7 (E) and T and also to 
get a better understanding of the underlying mechanism. It is worthwhile to mention that 
the predissociation lifetimes (r) are usually of the order of picoseconds and as a result there 
are only a few numerically exact TDQM studies available in the literature[167]. In our cal- 
culations we have employed the wave packet golden rule method[168] to determine T. This 
chapter is divided into three sections. In sections 2 and 3 we discuss the methodology and 
the results obtained for photo dissociation and predissociation processes. The conclusions 
are given in the final section. 

5.2 Photodissociation 

As mentioned in the Introduction, photoabsorption from the ground electronic state (X^II) 
to l^E", l^A, B^E"^ and 2^11 channels lead to direct photodissociation. Therefore, we have 
considered only these four states in our calculation. We have made use of the ab initio 
PE curves and the transition dipole moments reported earlier in the literature[282-284]. 
Potential energy values as a function of OH bond distance R, for the ground state, were 
fitted, using an extended Rydberg function[300]: 

V(r) = -De(l-0 + aiP + + azp^)exp{-')l p) + De (5-1) 

where p = R - Rg with Rg the equilibrium bond distance and De the classical bond 
dissociation energy. The parameters ai, a 2 , as and 7 obtained from the fit are given in 
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Table 5.1: Extended Rydberg function parameters for theX^R and A^I,+ potential-energy 
curves of OH. 



x^n 

A^E+ 

D,(eV) 

re(au) 
ai (au“^) 
a2(au"^) 
a3(au-3) 
7(au"’) 

4.47 

1.8519 

2.448 749 07 
1.609 992 94 
0.690 901 18 
2.481 589 59 

2.54 

1.9294 

3.026 579 60 
2.693 375 05 
1.120 864 54 
3.065 781 57 


Table 5.2: X^H state eigenvalues for u"=0-5. 


t/ 


E^(cm 


0 

1 

2 

3 

4 


1792.79 

5263.87 

8601.42 

11751.86 

14767.36 


5 


17597.10 
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table 5.1. The excited state PE curves and the transition dipole moments were interpolated 
smoothly using cubic splines. 

Once the ground electronic state PE curve is known, the next step is the computation 
of the ground vibrational eigenfunction. In order to compute this we have made use of 
the spectral method[83, 84] which consists of two steps. In the first step, vibrational 
eigenvalues corresponding to the groimd electronic state are computed and in the second 
step the vibrational eigenfunctions are computed. For the first step, we started with the 
initial wavefunction 

0) = .4 e~^ (5-2) 

where A — 0.4454, B — 0.0309 au ^ and = 3.5 au and time-evolved using the split 
operator method on the ground electronic state PE curve. At each step of the propagation 
the autocorrelation function 

c(t) = (V-(fi,0)|^(iJ,f)) 

was determined. We carried out the time evolution for 16,384 steps with each time step 
equal to 0.133 fs. After completing the time evolution, c{t) was Fourier-transformed to 
obtain the power spectrum, c{E) as shown in figure 5.2. For graphical resolution purpose 
we have actually plotted log c{E) in figure 5.2 which shows several sharp maxima and each 
one of them corresponds to a particular eigenvalue of the bound state. In order to obtain 
a better energy resolution, the time-evolution had to be carried out for a long duration as 
deternoined by 

AE = ^ (5-3) 

where T = NTx DTv/ith. NT and DT the total number of time-steps and the time-increment 
respectively. In addition, to avoid the overlapping sidelobes of the energy spectrum, we 
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multiplied the c{t) by Hanning window function w{t)/T, before Fourier transformation, 
as suggested by Feit et al[83]. Also the accuracy of the eigenvalues could be improved 
further by using a single line fitting procedure[83]. The eigenvalues obtained for the first 
six vibrational levels are listed in table 5.2. 

Once the eigenvalues were known, the vibrational eigenfunctions could be obtained 
from 

= N r dt (5.4) 

J-00 1 

where Ei is the eigenvalue of a particular vibrational level and N is the normalisation 
constant. 

Photodissociation cross section(cr) was computed as a function of energy E{ = Ei + 
hu, where u is the frequency of the incident radiation and Ei the internal energy of the 
molecule) using the formula[179, 224]: 

C(t) 

<t{E) = ^ r dt (5.5) 

3aC€o J-oo 

where eo is vacuiun permittivity constant. 

\<f>{t) >= e-’^“‘/''l<^(0) > (5.6) 

where 

l<;i(0) >= mIV’(O) > 

Here H,:, is the Hamiltonian of the excited electronic state and C{t) is the autocorrelation 
function defined as indicated in eq.(5.5). 

The initial wavefunction for the ground vib-electronic state, (i/»«»=o), obtained by the 
spectral method was transposed on each one of the excited electronic states after multiply- 
ing it with the corresponding transition dipole moment. Time-evolution of the tran p 
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wavefunction, ^(0), was accomplished by a second order differencing(SOD) scheme[85, 88] 
and the spatial derivative was obtained by the FFT method [85]. The spatial grid consisted 
of R = 1.5 (0.1) 14.2 au and we had to time evolve the wavefunction for states 1^E“, l^A 
and for a total of 8192 steps with each time step being, At = 0.00241889 fs. This 

means that the photodissociation through these states is over in about 8 fs. This can also 
be seen from the readily decaying autocorrelation function (|C'(t)|) for these states shown 
in figure 5.3a. 

The photo dissociation through the 2^11 state took significantly longer time; about 25 fs. 
We had to time evolve the wavefunction for 16,384 steps with At = 0.001584 fs. It is clear 
from figure 5.3a that \C{t) \ for the 2^11 state decays relatively slowly. An examination of the 
wavefunction reveals additional structures developing with time as illustrated in figure 5.4. 
As a matter of fact, the wavefunction spreads over the entire spatial grid and gets reflected 
from the edge. In order to avoid the reflections we used an absorbing potential [208]: 

y f ^ R ^ Robs ~f' ^Rabs f 8^ 

\ 0 ‘ R<R<j>s 

where Vjo is the hei g ht of the absorbing potential and ARabs its width. We have used V^o 
= 0.01 hartrees, R„6s = 10.5 au, ARabs = 3.5 au. The difference in time-evolution on the 
2^n state can be traced to the presence of a small hump in the PE curve which arises due 
to the strong interaction between 2^11 and 3^11 states (not shown), van Dishoeck et al[301] 
had suggested that this problem be handled using diabatic formulation but we have been 
able to solve it within the adiabatic representation. 

Values of cr(E) computed from the autocorrelation functions for the different excited 
states are presented in figure 5.3b along with the earlier theoretical [282-284] and exper- 
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Figure 5.3. (a) Autocorrelation function (|C(t)j) and (b) the photodissociation cross 
section (cr(E)) corresponding to l^E", l^A, and 2^11 states. 1 Mb = 1.0 x 10“^* cm^ 


14.5 
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Figure 5.4. Plots of |Vf as a function of R at different time steps for photodissociation 
j^om 2^n state. Each time step = 0.001584 fs. 
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imental [284, 285] results. While the maxima in o’(E), including an additional hump for 
the X^n-2^n excitation, are located at about the same place as in the time-independent 
quantal results, there is a noticeable discrepancy in the magnitude of the maviTnum a 
values. A careful examination of the results reported by van Dishoeck and Dalgarno[283] 
reveals that their results are dependent on the choice of basis set and hence differences of 
the kind noticed in figure 5.3 are not surprising. In any case, the fact that the two different 
theoretical approaches yield nearly the same cr(E) implies that the observed discrepancy 
between theory and experiment is likely to be due to some error in the experimental re- 
sults. One obvious source of error would be the uncertainty in the initial vib-rotational 
state of OH, in experiments. 

5.3 Predissociation 

It has been mentioned in the Introduction that the predissociation takes place in the first 
excited state, due to spin-orbit coupling between A^E"*" and ^E“, 1^E“ and “^H states 

(shown in figure 5.5). We have also mentioned that earlier investigations [287-298] on the 
OH predissociation suggested that the ^E“ state is responsible for the observed weak pre- 
dissociation near v’ < 2 and the ^H state is responsible for the strong predissociation near 
v’ > 5. Since the ^E“ and '^H states cross the A^E"^ state near v’ = 2 and 5 respectively, 
the above proposed mechanism is quite self-explanatory. However, when we want to com- 
pute r and its dependence on various vibrational levels of A^E"^ state, we have to consider 
not only the location of the crossing between the A^E"*" and the repulsive states but also 
the role of the spin-orbit coupling between the states and its i? dependence. Therefore, 
we have included the i?-dependence of the spin-orbit coupling explicitly in our present 
calculation. Since Lee and Freed [288] have already shown that the asymptotic coupling 
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between the repulsive states is not important while determining T, we have not considered 
the same in the present study. 

The ab initio PE curves and the spin-orbit coupling matrix elements reported by 
Yarkony[298] were connected smoothly through cubic spline interpolation. Since isolated 
resonance approximation has been shown to be valid for the present system[287], we have 
computed the contribution due to continuum of each repulsive state a, independently and 
the total r is obtained by summation, i.e., P*®* = 

As has been pointed out by Villarreal et al.[168], an exact time-dependent treatment of 
predissociation would involve inordinately long time evolution(see section 2.3.7). Therefore 
we followed the wave packet golden rule treatment suggested by them[168] for computing 

r^: 

r. = > * (5-9) 

where |(;6i(0) > = |^i(0) > and |<^f > = >• '0i(O) is the initial wave- 

function corresponding to the = 1,2, . - 7) vibrational level of the state. 

represents the spin-orbit coupling matrix element and Ha is the Hamiltonian correspond- 
ing to repulsive state. The time-evolution of the wavefunction on each repulsive state 
was accomplished by a second order differencing(SOD) scheme with the FFT routine[85] 
to obtain the spatial derivative of the wavefunction at any time t. The grid consisted of 
R = 1.25 (0.065) 9.505 au and the time-evolution involved a total of 8192 steps with each 
time step being 0.00126 fs. 

The ab initio PE values for the A^E''" state [298] were also fitted to an extended Rydberg 
function and the resulting parameters are included in table 5.1. The vibrational eigenvalues 
and eigenfunctions of the A^E"'' state were obtedned by the spectral method [82], using the 
computational scheme described in the earlier section for the X^H state. The vibrational 
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eigenvalues were obtained from the power spectrum shown in figure 5.6a, and they are 
given in table 5.3 for v’ = 0, 1, 2, . . . 7. The corresponding eigenfunctions are shown in 
figure 5.6b. We have verified that the spacing between vibrational levels of and X^II 
states computed by us are comparable to the available experimental results [294]. 

The computed values of ra(a = and ^11) are reported in table 5.4 and 

also plotted in figure 5.7 for different v’ along with the available experimental and time- 
independent quantal results. Our results reproduce qualitatively the trend in with 
increasing v’ reported earlier. That is, increases with increase in v’ until it reaches a 

maximum at v’=6 and declines subsequently. Although it is not apparent from figure 5.7, 
it is clear from table 5.4 that our estimate of T*®* for v’=l and 2 is substantially larger than 
that reported by others. The particularly low value of F obtained by Yarkony[298] can 
be explained readily as he considered only the predissociation into the state. Sink et 
al[287] had not included the variation of the spin-orbit coupling with R. The lower value 
of r obtained from experiments could be indicative of the reduced predissociation from 
the higher rotational states of OH. 

For v’=3, r*®* computed by us is compaurable to the value obtained from experiments 
and it falls exactly in between the two earlier theoretical estimates. For v’=4, our F*®* 
is close to the previous theoretical results, but is an underestimate when compared to 
experiment. For v’=5-7-, there are no experimental results available and our results are in 
excellent accord with the results of Sink et al[287]. 

It is worth emphasising that the results in table 5.4 demonstrate clearly that the pre- 
dissociation is dominated by the spin-orbit coupling with the state[116]. The increase 
in F‘°^ with increase in v’ through 6 can be easily understood from the location of the 
crossing points. As shown in figure 5.5, ^E“, l^E" and ‘^11 states cross the A^E"^ state 
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Table 5.3: state eigenvalues for t/=0-7. 


V 


E^(cin 


0 

1 

2 

3 

4 

5 

6 


1440.07 

4338.96 

7125.64 

9706.55 

12120.33 

14288.59 

16271.85 


7 


17974.65 


□ Sink etal.[287] 

O Copeland etal.[295] 
X Czarny etal. [291] 


Figure 5,7. Predissociation linewidths (Pa) for each channel and the T as a function of 
v’. Some of the earlier results are indistinguishable from the ones shown for v’=l-3 and 
therefore are not shown separately. 
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Table 5.4: Predissociation linewidths,r (in cm”^) computed using the golden rule treat- 
ment. 

The numbers in parentheses represent the powers of 10. 


v’ 



“n 

"ptcft 

r(others) 

Ref. 

1 

8.2(-6) 

3.6(-6) 

7.6(-4) 

7.7H) 

9.0(-8) 

[6a] 






6.6(-13) 

[298] 






< 2.1 (-7) 

[7b] 

2 

2.9(-5) 

7.7(-6) 

8.4(-4) 

8.7H) 

2.3(-5) 

[6a] 






3.3(-5) 

[298] 






1.8(-5) 

[7a] 






2.9(-5) 

[7b] 

3 

1.9(-2) 

4.1 (-5) 

9.6(-4) 

2.0(-2) 

1.2(-2) 

[6a] 






2.4(-2) 

[298] 






1.9(-2) 

[7d] 






1.8(-2) 

[7e] 

4 

1.7{-1) 

1.5(-2) 

2.9(-2) 

2.1(-1) 

1.9(-1) 

[6a] 






1.9(-1) 

[298] 






«7.5(-l) 

[7c] 

5 

5.4(-2) 

1.4(-1) 

8.3(-l) 

1.0 

1.0 

[6a] 






2.1(-1) 

[291] 

6 

3.9(-2) 

4.4(-2) 

1.5 

1.6 

1.8 

[6a] 






1.5 

[291] 

7 

4.9(-2) 

5.6(-2) 

3.1(-1) 

4.2(-l) 

4.5(-l) 

[6a] 






1.5(-1) 

[291] 
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near the vibrational levels v’=3-6 corresponding to J? = 2-3 au for which increases 
with R. For iZ > 3 au, the for the dominant channel, ^11, decreases with increase 
in R explaining the decline in beyond v’=6. For the and the ^E“ channels 
continues to increase with increase in R thus explaining the continued increase in F in 
going from v’=6 to v’=7. 

5.4 Conclusion 

The TDQM approach has been applied to study both photodissociation and predissociation 
processes in OH. The photodissociation cross section values computed by us are comparable 
to the time-independent quantal results but a factor of two smaller than the experimental 
results. 

The predissociation linewidths calculated using the Fermi golden rule treatment are in 
general accord with the earlier theoretical results. In general, it has been shown that F*®* 
increases with increase in v’ through 6 and then it declines subsequently and that this is 
explainable on the basis of the R-dependence of the for the donainant channel ^H. 



Chapter 6 


Summary and Conclusion 

Time-dependent quantiun mechanical approach is increasingly being used in recent 
years to investigate various elementary chemical processes such as photodissociation, pre- 
dissociation, reactive scattering, molecule-surface scattering etc.. In the present thesis the 
TDQM approach has been employed to study (i) photo-predissociation of Nal, (ii) pho- 
todissociation of HI and (iii) photodissociation and predissociation in OH radical. Impor- 
tant conclusions derived from oiu studies are summarised below. 

In chapterl of the thesis, we have reviewed the necessary background material related 
to the TDQM approach to photodissociation and predissociation processes. We have also 
presented a short literature review emphasising mainly on the recent studies. 

Theoretical and computational methodologies related to the TDQM approach to pho- 
todissociation and predissociation processes have been presented in chapter 2. This chapter 
has been divided into three sections. In section 1, we have described the existing compu- 
tational methodologies for implementing the TDQM approach. Section 2 of this chapter 
is devoted mainly for the interaction of light with molecules. We have presented a semi- 
classical and a quantum mechanical description of the photoabsorption process. We have 
also described, how the photodissociation cross section is computed in the time-domain. 
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Also we have presented a critical analysis of the nature of the wavefunction prepared by 
a laser pulse on the excited electronic state. In section 3, we have described the computa- 
tional methodologies for investigating multieletronic state problems. Also the golden rule 
treatment for computing the predissociation linewidth is included in this chapter. 

In chapter 3 of the thesis we have presented the results of the TDQM study of the 
photo-predissociation of Nal using ultrashort laser pulse. We have investigated the nature 
of the wave packet prepared by a femtosecond laser pulse, whose excitation frequency is 
chosen near, very near and far from resonance excitation frequency. Since the excita- 
tion frequencies are different, the WPs prepared on the excited electronic state also have 
different characteristics. The norm of the WP on the excited state during the pulse prepa- 
ration step has been analysed. After the pulse preparation step is over, the WP is time 
evolved for about 3 ps on the covalent excited electronic state which is coupled to the 
ground electronic state. Due to the diabatic coupling between the ionic ground state and 
the covalent excited state, the time evolving WP bifurcates near the crossing point and a 
s mall fraction of it (as determined by the Landau-Zener curve crossing probability) leaks 
through the crossing to form neutral products Na and I and the remaining part of the 
WP oscillates between ionic and covalent turning points. Whenever the WP reaches near 
the crossing point a small fraction of it leaks through crossing. We have estimated the 
oscillation periods corresponding to photoexcitation wavelengths 300, 321 and 344 nm and 
our results are in good agreement with the experimental results of Zewail and coworkers. 
However, when the excitation wavelengths aie far from resonance, i.e., Aex = 360 and 390 
nm, the WP spreads immediately ( in less than 0.5 ps ) and shows several interference 
pattern. Therefore, oscillation periods estimated by us are not in agreement with experi- 
mental results. We have observed that the transient effect and the energy of the WP are 
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highly responsible for the unusual characteristics of the WP evolution. 

In chapter 4 of the thesis, a detailed analysis of the photodissociation of HI in its first 
absorption band has been made. HI photodissociation takes place around 250nm resulting 
in H + I and H + I* formation. In this energy range, the following potential energy curves 
participate in the photodissociation: ^Hi, ^Hi, and ^Hq. Among these, the former two 
correlate with the H + 1 channel and the later two correlate with the H + 1* channel. Except 
for the ^Ho state, all other states are non-adiabatically coupled to each other. We have 
computed the photodissociation cross section(<7) corresponding to the photoexcitation firom 
the ground vibrational level of the ground electronic state to all of the above mentioned 
excited states. We have made use of the PE curves, non-adiabatic coupling potentials and 
the coordinate dependent transition dipole moments reported earlier by Levy and Shapiro. 
Our computed a and P values are in excellent agreement with the experimental and the 
time-independent quantum mechanical results. We have also computed the anisotropy 
parameter(,5) for the H -t- P channel and our values as a function of excitation energy 
are in accord with the experimental and the time-independent quantal results. In addition, 
we show that the P value changes significantly with vibrational excitation of the parent 
molecule and point out vibrational excitation prior to photodissociation as a means of 
channel control. The results are interpreted in terms of the influence of the non-adiabatic 
coupling on the time-evolution of the wavefunction. 

In chapter 5, we have described our results obtained for the photo dissociation and 
predissociation processes in OH radical using the TDQM method. Since OH radical is 
one of the important constituents in the astrophysical environment, understanding various 
dissociation pathways in OH becomes particularly essential. Photodissociation in OH 
takes place due to electronic transition from the groimd electronic state (X^H) to the 
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excited electronic states, l^E", l^A, and 2^11. The PE curves and the coordinate 

dependent transition dipole moments are taken from the ab initio results of Dalgarno and 
coworkers. Our computed photodissociation cross sections are in good agreement with the 
time-independent quantal results and a factor of two less than the experimental results. 
Predissociation in OH takes place from its first absorption band (X^n-A^E"^) due to spin- 
orbit coupling between the bound state A^E'*' and the repulsive states ^E“, 1^E“ and '*11. 
Using the recent ab initio results of PE curves and the spin-orbit matrix elements reported 
by Yarknoy we have obtained predissociation linewidths corresponding to predissociation 
from v’=l-7 of the A^E"*" state. The results are analysed in terms of the crossings between 


the PE curves. 
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